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ABSTRACT 



The Navy Operational Global Atmospheric Prediction System 
(NOGAPS) boundary layer and cloud parameterizations are 
evaluated for a case of explosive cyclogenesis in the 
western North Pacific Ocean. Storm-region diabatic heating 
estimates are obtained from the application of quasi- 
Lagrangian diagnostics to a thermodynamic energy budget 
calculation, and from the NOGAPS diagnoses. 

The model-diagnosed sensible heating appears to be 
correctly positioned, while the diagnoses of convective, 
large-scale and open-cell cumulus condensation heating 
produce cloud features which generally reflect the distri- 
bution of the clouds in the satellite imagery. 

NOGAPS provides a better estimate of the diabatic heat- 
ing over the open ocean than does the thermodynamic energy 
budget calculation. The contribution of diabatic heating in 
this case study is determined to be as significant as that 
of thermal advection. This diabatic energy input is an order 
of magnitude greater than that determined from studies of 
continental cyclogenesis. 
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INTRODUCTION 



I . 

Explosive extratropical maritime cyclogenesis poses 
extreme hazards to naval and commercial shipping as well as 
to stationary oceanic platforms. Rapid cyclogenesis is 
characterized by intense deepening of a storm system over a 
short period of time, which results in the generation of 
high wind and damaging waves. Present numerical models 
often fail to forecast the intense nature of these systems, 
which leads to storm-related property damage, injury and 
death. 

The explosively deepening cyclone is described by Sanders 
and Gyakum (1980) as primarily a wintertime marine event. An 
explosive event is defined when the normalized (sin (p/sin 60) 
central pressure fall is equivalent to one mb/h for 24 h at 
60 °N. Their study indicates that these cyclones occur in 
the vicinity of a strong sea-surface temperature gradient 
along the leading edge of an outbreak of cold, continental 
air. The regions of maximum frequency were found to be in 
the western North Pacific and Atlantic Oceans. 

The failure of numerical models to satisfactorily fore- 
cast the rapid intensification of explosive systems is well 
documented in several studies, including those of Sanders and 
Gyakum (1980) , Bosart (1981) , and Gyakum (1983) . These 
studies elucidate the problems of numerically forecasting 
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explosive cyclogenesis. Sanders and Gyakum verified the 
National Meteorological Center (NMC) six-and seven-layer 
primitive equation (PE) model predictions of explosive 
events observed during the 1977-1979 seasons. Their results 
demonstrated that the coarse mesh, six-layer PE model forecast 
only about 25% of the observed 12-h central pressure tendency 
during the explosive stage of cyclogenesis. Use of the higher 
resolution seven-layer PE model led to only a slight improve- 
ment in the forecast (< 10%) . Bosart (1981) documented simi- 
lar deficiencies with the Limited-area Fine Mesh II (LFMII) 
model in his study of the 1979 President's Day storm. Again, 
the magnitude of the central pressure fall was drastically 
under forecast , and this intense storm struck the eastern 
seaboard of the U.S. with little warning. 

In an extreme case, Gyakum (1983) found that the LFMII 
model underforecast the 12-h explosive deepening of the 1978 
Queen Elizabeth II storm by 55 mb. This failure led to gross 
errors in the wind and wave forecasts and resulted in damage 
to commercial shipping, most notably in the Queen Elizabeth 
II, for which the storm was named. Calland (1983) also dis- 
covered significant forecast errors in both storm intensity 
and track in his evaluation of the NMC and Fleet Numerical 
Oceanography Center (FNOC) coarse-mesh primitive equation 
model predictions for a western North Pacific Ocean system. 

In numerical simulations of the Queen Elizabeth II storm, 
Anthes et al. (1983) found the incipient model cyclone to 
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be shallow and sensitive to the initial conditions. They 
suggested that increased vertical resolution in the lower 
troposphere and an initial analysis with more representative 
moisture, static stability and low-level wind structures 
will lead to improved operational forecasting of maritime 
cyclogenesis . 

This thesis is part of a larger investigation into the 
nature of maritime extratropical cyclones, which has an 
overall objective of improving numerical weather prediction 
over the oceans. The primary objective of this study is to 
evaluate the contribution of diabatic heating to storm 
development, especially during the explosive deepening stage. 
A secondary objective is to analyze the performance of the 
cloud and boundary layer parameterization package of the UCLA 
General Circulation Model as incorporated in the Naval Opera- 
tional Global Atmospheric Prediction System (NOGAPS) . 

Following the work of Calland (19 83) , an explosive North 
Pacific Ocean cyclone is analyzed to determine the quantita- 
tive contribution of diabatic heating to cyclone development. 
This storm possesses the characteristics of 'Type A' cyclo- 
genesis which is described by Petterssen et al . (1962) and 

Petterssen and Smebye (1971) as frontal wave development 
initiated by baroclinic instability. This type of cyclo- 
genesis usually occurs with nearly straight upper-level 
flow over a zone of maximum baroclinity. As the storm 
develops, an upper-level cold trough forms and maintains a 
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nearly constant spatial relationship with the low-level cy- 
clone. Petterssen indicates that thermal advection is a 
dominant term in 'Type A' cyclogenesis, while vorticity 
advection plays only a minor role. A detailed discussion of 
the synoptic development of the storm analyzed in this the- 
sis is given by Calland (19 83) . The significant events of 
that development are summarized in Appendix A to present a 
profile of cyclogenesis as it occurred between 13-15 January 
1979 . 

The quantitative analysis of the effect of diabatic heat- 
ing on storm development is accomplished using the NOGAPS 
heating package (C0MP3) in conjunction with the quasi- 
Lagrangian diagnostics routine as developed by Johnson and 
Downey (1975 a and b) and applied by Calland (1983) in his 
mass and circulation budget study. Fields of convective 
condensation heat release, large-scale condensation heat 
release, sensible heat flux, long wave radiation heating, 
short wave radiation heating and total diabatic heating are 
derived from diagnostic application of the NOGAPS heating 
package to the European Center for Medium-Range Weather 
Forecasting (ECMWF) level Illb analyses of the First Garp 
Global Experiment (FGGE) data base. This diagnostic use of 
NOGAPS is similar to the semi-prognostic model application 
described by Lord (1982) , which is discussed in the next 
chapter. The FGGE data base and ECMWF analyses are discussed 
in Appendix B . 
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The ECMWF analyses are interpolated to the six sigma 
levels of the model at 12 h intervals for a grid which 
covers the region of storm development with a horizontal 
resolution of 1.875° in latitude and longitude. These 
interpolated analyses are used with the NOGAPS diabatic 
package to diagnostically estimate the heating fields. 

Storm volume averages of each heating field are prepared for 
the immediate area of storm influence. Diabatic energy 
inputs to the column from the total heating and latent heating 
fields are estimated for the storm volume. 

In a parallel analysis, the FGGE data base is used to 
estimate the total and advective temperature tendencies in 
the local storm environment through the application of 
quasi-Lagrangian diagnostics. In this application, the storm 
volume is centered on, and translates with, the cyclone. 

This diagnostic technique effectively removes developmental 
aspects of the cyclone which are associated with the storm 
motion. Another estimate of the diabatic temperature change 
is derived as a residual between the total and advective 
temperature change fields. These results are compared with 
those obtained from NOGAPS as a check on the model's 
effectiveness . 

A survey of the literature on the role of surface fluxes 
and latent heat release in cyclogenesis is included in 
Chapter II. The boundary layer and cloud parameterizations 
are evaluated in Chapter III through diagnostic use of the 
NOGAPS code. S torm-rel ated diabatic heating, as determined 
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by the application of quas i-Langrangian diagnostics to the 
FGGE data, is compared to the model-generated diabatic 
heating in Chapter IV. Conclusions and recommendations for 
future study are incorporated in Chapter V. 



17 



II. LITERATURE SURVEY 



A. GENERAL 

The surface fluxes from the ocean and the release of 
latent heat through cumulus convection play important roles 
in ocean cyclone development. When these processes occur in 
the proper location with respect to a developing system, 
they can enhance cyclone growth. A review of past studies 
suggests that these diabatic processes are important 
ingredients in explosively developing extratropical maritime 
cyclones . 

B. SENSIBLE AND LATENT HEAT FLUX 

Petterssen et al. (1962) constructed composite models of 

North Atlantic cyclones at various stages of development. 

They documented strong fluxes of sensible and latent heat in 

the cold, southward flowing polar air to the north and west 

of a well-developed frontal zone over the ocean. Although 

both fluxes were strong in this region, sensible heating was 

dominant. They found that the daily loss of heat from the 

ocean to the atmosphere in extreme cases could exceed 20 

times the average absorption of incoming solar radiation at 

high latitudes during winter. The typical sensible heat 

2 

flux maximum was approximately 700 W/m . The typical latent 
heat flux maximum was found to be of this same order of 
magnitude. Petterssen indicated that the addition of the 
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sensible heat and moisture flux to cold, polar air will 
destabilize the boundary layer and cause some convective 
activity to the north of the frontal zone. Although the 
sensible heating may be extreme in this region, deep convec- 
tion does not occur unless the air curves cyclonically . 

These composite models depict the warmer air to the south 
of the frontal zone as tropical in nature. South of the 
frontal zone the sensible heat flux is minimal, since the 
atmosphere is essentially in thermal equilibrium with the 
ocean surface. Most of the moisture that is provided to the 
frontal system, and later removed as precipitation, has 
tropical origins . 

Anthes et al. (1983) determined through model simulations 
of the Queen Elizabeth II storm that the surface fluxes have 
a moderate effect on storm development. In agreement with 
Petterssen, they found large changes in the PBL and low- 
level frontal structure to the southwest of the cyclone in 
the area of strong cold advection. The simulations indicated 
that the addition of the surface fluxes reduced the gradient 
of the front as the cold air was modified by the warmer ocean 
At the same time these fluxes destabilized and moistened the 
boundary layer, especially to the southwest of the storm 
center . 

Based on a sample of 256 events of explosive cyclogenesis 
Sanders and Gyakum (1980) placed the location of 92% of these 
storms within, or to the north of the maximum baroclinity 
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and wind at 500 mb. This statistic indicates that most of 
the explosive events were placed in the polar air mass 
where the large-scale horizontal temperature contrasts and 
the transfer of sensible and latent heat from the warmer ocean 
to the colder air are important. Most of these explosive 
deepeners occurred during winter with only a few summer cases 
noted. This seasonal dependence supports the conclusion that 
large-scale, horizontal temperature contrasts are important 
to explosive development. Mullen (1983) also found a strong 
preference for explosive deepening to the north of the maxi- 
mum baroclinity at 500 mb in his investigation of explosive 
cyclones in polar airstreams. 

Sanders and Gyakum found that extratropical explosive 
cyclones are not as sensitive to the sea-surface temperature 
as tropical cyclones. Explosive cyclones occur over a large 
range of sea-surface temperatures (0°-23°C) . The sensible 
and latent heat exchange at the surface is most intense when 
fast-moving, cold air crosses a strong sea-surface tempera- 
ture gradient toward warmer water. This scenario results 
in a modification of the air mass and produces lower static 
stabilities in the lower troposphere, which is an important 
factor in explosive development. 

Staley and Gall (1977) demonstrated that growth rates of 
short baroclinic waves were enhanced by lower static sta- 
bility and strong low-level vertical wind shear. Sanders 
and Gyakum (1980) hypothesized that the positions of strong 
sea-surface temperature gradients would indicate areas of 
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low-level baroclinity with a higher susceptibility to explo- 
sive development. Indeed, their analysis of explosive 
developments in both the North Pacific and Atlantic Oceans 
demonstrated that these cyclones tended to develop in regions 
of strong sea-surface temperature gradients. 

Bosart (1981) also emphasized the important role played 
by surface fluxes of sensible and latent heat in the vicinity 
of a strong sea-surface temperature gradient. In his analy- 
sis of the President's Day storm, Bosart noted that large 
sensible and latent heat fluxes in the region of the Gulf 
Stream strengthened a low-level baroclinic zone and caused 
conditionally unstable lapse rates. In a later study of 
this same system, Bosart and Lin (1984) calculated the sensi- 
ble heat flux and found that the maximum ranged between 

2 

approximately 400 and 800 W/m during the period of cyclone • 

development. Latent heat fluxes were found to be as much 

as three times greater. The average of sensible heating 

estimated by Chou and Atlas (1982) during the initial cold 

air outbreak associated with the President's Day storm was 
2 

249 W/m . These large boundary layer fluxes altered the 
environment and enhanced the potential for deep convection, 
which provided the cyclone circulation with the heat and 
moisture required for explosive deepening. 

In a North Atlantic cyclogenesis case study by Gall and 
Johnson (1971) , the surface sensible heat fluxes generated 
significant amounts of available potential energy during the 
period of cyclone development. They indicated that forcing 
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of large-scale cyclogenesis was influenced by the horizontal 
extent of the transfer of energy at the surface over the 
ocean, which was determined by the large-scale environment. 

They found that the total diabatic available potential 
energy generation within the storm scale was sufficient to 
offset the frictional dissipation of kinetic energy during 
the early and mature stages of the cyclone. This comparison 
illustrates the important role played by sensible heating, 
as well as the other diabatic heating sources, in cyclogenesis. 

The direct dynamic contribution of sensible and latent 
heat flux to explosive cyclogenesis appears to be small. 

Gyakum (1983) cites the importance of low-level baroclinic 
forcing in initiating cyclone development and providing the 
lifting necessary to establish a CISK-like mechanism which 
could drive explosive growth rates. The direct contribution . 
of the surface fluxes in a dynamic situation of this type 
appears to be minimal, but Gyakum feels that these fluxes 
are important in that they destabilize the boundary layer 
and provide the air converging into the low center with the 
moisture necessary to sustain deep convection. 

C. LATENT HEAT RELEASE 

The role of latent heat release in cyclone development 
was first studied early in the nineteenth century. Margules 
(1903) calculated that in a potentially stable situation the 
release of latent heat served mainly to decrease the rate of 
cooling of a rising air parcel but did not act to increase 
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the parcel's kinetic energy. This study did not take into 
account the effect of the surrounding environment on the 
system or the rate at which the energy conversion took place. 
The conclusions of Margules ' study did, however, influence 
research in cyclogenesis for many years (Danard, 1964) . 

It was not until the 1950's that the convective process 
was again viewed as an important mechanism in cyclone develop- 
ment. Aubert (1957) concluded that released latent heat 
greatly increases vertical motion and tends to lower pressure 
surfaces in the lower troposphere and raise them in the upper 
troposphere. Petterssen et al . (1962) described the role of 

released latent heat in cyclone development as having two 
separate effects depending on the location of the released 
heat relative to a frontal zone. If the latent heat release 
• takes place in the cold, polar air to the north of a well- 
developed frontal zone, it has an effect much like that of the 
sensible heat flux. Condensation here serves to reduce the 
baroclinity as well as the potential energy of the cyclone 
as compared to a purely adiabatic system. If, however, the 
release of latent heat occurs over the frontal surface, 
condensation serves to increase the potential energy of the 
cyclone . 

Danard (1964) calculated the effect of released latent 
heat on vertical velocity when the static stability is 
allowed to vary horizontally. This allowance was an improve- 
ment over Margules' closed boundary condition. The results 
demonstrated that the effect of released latent heat on a 
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rising air parcel is to increase its vertical velocity. 

Upward motion is amplified in regions of heavy precipitation, 
with slight downward motion occurring in the surrounding 
areas. The increase in vertical velocity in the mid- 
troposphere from latent heat release leads to increased 
convergence at low levels and increased divergence aloft. 

Danard also related this release of latent heat to the 
production of kinetic energy and vorticity at upper and 
lower tropospheric levels which is of the same order of 
magnitude as the energy and vorticity production from dry 
adiabatic processes. Danard described the role of released 
latent heat as a mechanism for cyclone intensification, 
since he felt that a pre-existing disturbance had to be 
present for the necessary deep convection to occur. This 
description is supported by Mak (1982), who found that baro- 
clinic forcing organizes latent heat release on a scale 
similar to that of the existing disturbance. 

Danard (1966) modeled the vertical distribution of latent 
heat release as a parabolic function with a maximum in the 
middle troposphere. He compared the results of an energy 
analysis case study to calculated changes in potential energy 
due to latent heat release. He found that the increase in 
available potential energy from condensation heating was 
much more significant than the loss caused by the increased 
vertical motion in the low and middle troposphere. This re- 
sult led to the conclusion that release of latent heat would 
reduce the loss of available potential energy of the system 
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normally associated with vertical circulations. He calculated 

2 

a storm area latent heat release energy maximum of 11 W/m . 

Bullock and Johnson (1971) used a translating budget 
volume to estimate the generation of available potential 
energy for early, mature and occluding stages of an extra- 
tropical cyclone over the continental United States. They 
found that the horizontal and vertical distribution of the 
latent heat release is important in determining its contribu- 
tion to storm generation. They determined that maximum heat- 
ing in the lower levels, especially when it occurs over the 
warm sector, produces the maximum generation of available 

potential energy. The estimate for storm volume latent heat- 

91 2 

ing* potential energy production was 1.31 x 10“ ergs/s (8 W/m ) , 

which is consistent with Danard's findings. 

Tracton (1973) evaluated the role of cumulus convection 

in cyclone development using numerical model forecasts. He 

tested the following hypothesis: 

In some instances of extratropical cyclogenesis, 
cumulus convection plays a crucial role in the initia- 
tion of development through the release of latent heat 
in the vicinity of the cyclone center. In such cases, 
dynamical models that do not adequately simulate 
convective precipitation, especially as it might occur 
in an environment that is unsaturated, will fail to 
properly forecast the onset of development. 

His evaluation of the NMC six-layer primitive equation 

model, NMC limited-area, fine-mesh model and FNOC five- layer 

primitive equation model forecast performance in 19 separate 

cases of continental United States cyclogenesis supported 

the hypothesis. In all cases where convection occurred near 
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the low center in an unsaturated environment, the numerical 
models failed to properly forecast cyclone development. 

In contrast to Danard, Tracton viewed convection as an 
initiation process in cyclogenesis when it occurred near the 
surface low center. He suggested that released latent heat 
caused these cyclones to develop earlier than they would 
have with only large-scale baroclinic processes affecting 
them. This inference was demonstrated by associating the 
model forecast lag time with cyclones where the vorticity 
advection was weak. He found that the largest deepening at 
the storm center occurs when the ratio of upper tropospheric 
to lower tropospheric latent heating is smallest. 

In a study of east coast cyclogenesis, Danard and 
Ellenton (1980) supported the conclusion of Bullock and 
Johnson (1971) . Danard and Ellenton found that it is the 
configuration of the latent heat release that determines 
this diabatic influence on cyclone intensification. They 
discovered that the Laplacian of the surface heating did not 
lead to intensification over the low center, even though 
surface fluxes were strong in the cold air behind the storm. 
In agreement with Gyakum (1983) , Danard and Ellenton felt 
that the surface fluxes contributed to the development of 
favorable vertical distributions of temperature and moisture 
which enhance later cyclone intensification. 

Sanders and Gyakum (1980) also studied the performance 
of the NMC primitive equation models in explosive situations. 
They found that the primitive equation models drastically 
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underforecast the intensity of developing cyclones. An 
increase in horizontal resolution of the model did not cor- 
rect the magnitude of the observed deficiency. Sanders and 
Gyakum concluded that several factors are poorly represented 
in current models. Based on the extent of cumulus convection 
observed in satellite imagery of explosive cyclones, they 
indicated that a better representation of the effect of cumulus 
convection and latent heat release, as well as an improved 
PBL formulation, is required for the model calculations. 

Gyakum (1983) found that the Queen Elizabeth II cyclo- 
genesis was initiated by shallow baroclinic forcing in a 
region of widespread potential instability, which led to the 
outbreak of cumulus convection in the vicinity of the sur- 
face low center. He observed that adiabatic, quasi-geostrophic 
dynamics could not account for the observed intensity of cy- 
clone development during its explosive stage. 

The following relationship was used by Gyakum to deter- 
mine the contribution of diabatic heating to storm development: 



3h 

3t 



C ♦ Vh 



V • Vh + 
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dt ' adiab ^3t ; diab 
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(1) (2) (3) (4) (5) 



where (1) is the mean layer thickness change following the 
storm center; (2) is the thickness change produced by the 
cyclone's movement; (3) is the change due to horizontal 
temperature advection; (4) is the change from adiabatic 
warming and cooling due to storm-scale vertical motion; and 
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(5) is the diabatic contribution to thickness change. The 
contribution of diabatic heating to the thickness change 
over the low center was determined through a residual calcu- 
lation using the parameters of observed thickness change, 
temperature and assumed vertical velocity profiles. The 
residual warming of the column over the surface low center 
was found to be very significant (about 16°C) during the 
12-h explosive period of this cyclone. 

Gyakum concluded that much of the converging air at the 
surface ascended within convective towers near the low. This 
ascending air provided the bulk heating that likely forced 
the vertical motion not accounted for by adiabatic, quasi- 
geostrophic dynamics. The vertical heating (q) profile was 
modeled with a mid-level maximum and zero values at the 
surface and top of the atmosphere. Linear profiles were 
assumed between these points. By varying the level of maxi- 
mum heating, Gyakum found that heating profiles having a 
maximum at low levels could force surface convergence and 
vertical ascent on a scale consistent with the observed 
extreme geopotential height falls. He calculated that 
cumulus- induced subsidence warming in the vicinity of the 
low during the explosive stage produced column warming of 
16 °C over a 12-h period. Warming of this magnitude was 
sufficient to account for the observed 12-h thickness changes 
in the column. 

Model simulations of the Queen Elizabeth II storm by 
Anthes et al . (1983) support most of Gyakum ' s findings. 
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Comparisons of simulations with and without latent heating 
produced upper tropospheric temperature differences of 8.9°C, 
with the warmer values associated with the latent heat re- 
lease from deep convection. The simulation with latent 
heating included was colder in the lower troposphere where 
the heaviest precipitation occurred, but warmer in the sur- 
rounding areas due to the compensating subsidence. The 
genesis in the model was not significantly affected by 
latent heat release during the early stages. Anthes , et 
al . (1983) indicated that the relatively weak contribution 

of latent heating during this early period supported the 
hypothesis that baroclinic instability was the mechanism 
for early development. Once stronger vertical motions were 
established, latent heating played a more important role in 
later intensification. 

In a study of an extratropical maritime cyclone during 
the initial phase of the 1975 Air Mass Transformation 
Experiment, Chen et al. (1983) also determined that maritime 
cyclogenesis is enhanced by latent heat release. Model 
simulations of this storm indicated that the latent heating 
was linked to the surface fluxes. When the surface fluxes 
were removed from the experiment, the amount of latent heat- 
ing, as well as cyclogenesis, was reduced. The latent heat- 
ing also had an impact on the phase speed of the cyclone. 
Reductions in latent heating led to a slower translation 
of the storm system. 
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In summary, numerous studies have found surface fluxes of 
sensible and latent heat, as well as the effects of latent 
heat release induced by convective processes, to be signifi- 
cant in certain situations of cyclogenesis. The surface 
fluxes provide heat and moisture to the cold continental 
air which destabilizes the boundary layer and provides a 
favorable environment for sustained convection. Latent heat 
release, particularly in the vicinity of a surface low 
center, will increase both the vertical velocity and surface 
convergence in the region. The lower the level where the 
maximum of this heat release occurs, the greater the surface 
convergence. This type of bulk heating will lead to an in- 
crease of the energy and vorticity of a cyclone, especially 
if the deep convection takes place near a surface low center. 
The latent heat release appears to be a mechanism for storm 
intensification rather than initiation, and increases in 
importance once the vertical circulation has been established 
by low-level baroclinic forcing. 

Several studies suggest that a CISK-like mechanism may 
be responsible for the rapid growth and extreme deepening 
observed in explosively deepened cyclones. Subsidence warm- 
ing over the surface low would be induced by surrounding 
convective towers and may provide enough heating to the column 
to produce the 12-h geopotential height falls observed in 
some of the intense warm-core systems. 

Numerical models that do not properly consider the 
diabatic effects in cyclone development perform poorly where 
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convection occurs in an unsaturated environment. All of the 



operational models that have been tested fail to forecast 
adequately the rate of development of explosive maritime 
cyclones . 
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III. NOGAPS DIABATIC PACKAGE EVALUATION 



A. GENERAL 

The Naval Operational Global Atmospheric Prediction Sys- 
tem (NOGAPS) model, which is based on the UCLA General 
Circulation Model (GCM) , is comprised of a series of routines 
which are structured to produce analyses and predictions on a 
global basis. A discussion of NOGAPS (version 2.0) is con- 
tained in Appendix C. Subroutine C0MP3 incorporates the 
main diabatic processes for the model. 

Diabatic processes of this Pacific storm system are 
analyzed using COMP3 in a diagnostic mode. This approach is 
similar to the semi-prognostic technique utilized by Lord 
(1982) in his test of the Arakawa-Schubert cumulus parameteri- 
zation for the tropical Atlantic region. Lord made use of 
observed parameters to estimate the large-scale forcing of a 
cumulus ensemble at given observation times. These observa- 
tions and the parameterization are combined to predict 
precipitation and cumulus warming and drying for each obser- 
vation time. The semi-prognostic method is not an integration 
in time and it is therefore free of any modeling errors other 
than those in the parameterization. This approach allows 
for a comparison of model-generated tendencies with observed 
tendencies for each observation time, since the observed ten- 
dencies are not included in the model calculations . The 
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diagnostic use of C0MP3 is described in more detail in the 
next section. 

Diagnostic (semi-prognostic) fields of sensible heat 
flux, convective condensation heating, large-scale condensa- 
tion heating, total surface precipitation, boundary layer (BL) 
thickness and layer cloud instability (LCI) are produced by 
C0MP3 using the ECMWF FGGE data as the large-scale forcing. 

The diabatic processes, as depicted by COMP 3 using the semi- 
prognostic technique, are evaluated during storm development 
by comparison with observed data and satellite imagery. 

B. PROCEDURE 

The 12-h FGGE data are utilized to study storm diabatic 
processes for a case of Western Pacific explosive cyclo- 
genesis during the period 00 GMT 13 January to 00 GMT 15 
January 1979. A complete discussion of the synoptic develop- 
ment of this cyclone is given by Calland (1983) . A summary 
of this development is included as Appendix A. 

COMP3 is evaluated independently of the full NOGAPS model. 
Boundary layer (BL) information which is normally provided to 
COMP3 by NOGAPS is developed internally through an iterative 
process at each 12-h evaluation time. Diagnostic (semi- 
prognostic) utilization of C0MP3 is accomplished by resetting 
the basic parameters of surface pressure, temperature, mois- 
ture and wind components to the ECMWF analysis values after 
each iteration. This technique allows COMP3 to compute the 
physical processes without advancing in time. Ten iterations 
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of the code are performed at each analysis time to allow the 
opportunity to develop a suitable boundary layer. 

An initial guess of 100 mb is made for BL thickness. The 
initial jumps of the wind components, temperature and mois- 
ture at the BL inversion are set equal to zero at 00 GMT 12 
January 1979. The iterated boundary layer parameters are 
passed to the next analysis time as the initial guess and 
the code is again iterated to develop the boundary layer for 
each new time period. This process is repeated for each time 
through 00 GMT 15 January 1979. It should be noted that a 
shift in the FGGE data- grid occurs at 12 GMT 14 January 1979 
(Calland, 1983) . For continuity, BL parameters in the over- 
lapping regions of the two data grids are passed directly 
to the new grid as the first guess field for 12 GMT 14 
January. Initial BL parameters for the non-overlapping ar'ea 
of the new grid are given by the northern row and western 
column of the overlapping portion of the grid. Due to 
oscillatory behavior in some of the BL parameters, averaging 
is performed over the last six iterations of each evaluation. 
More details on this boundary layer treatment are presented 
in the next section. 

Sea-surface temperatures (SST's) from FNOC are used from 
12 and 16 January 1979. The 12 January 1979 analysis corres- 
ponds to the first FGGE data grid and is used for the first 
three evaluation times. The 16 January 1979 SST field corres- 
ponds to the second FGGE data grid and is incorporated into 
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the analysis of the last two evaluation times. These analyses 
and the storm track are included as Fig. 1. The SSTs are 
extracted from the FNOC analysis at the FGGE gridpoints. 

These values are used as surface temperature in C0MP3, except 
over land where the FGGE 1000 MB air temperature is entered 
as the surface temperature. All fields discussed in this 
chapter utilize a 53 x 30 section of the ECMWF analysis grid 
with 1.875° latitude and longitude resolution. 

The diagnostic output is analyzed in the following sec- 
tions to demonstrate model behavior and performance in a 
case of maritime explosive cyclogenesis. Comparisons with 
the FGGE analyses, satellite imagery and previous investi- 
gations are made to determine the realism of the model's 
diabatic parameterization response in this intense developing 
cyclone . 

C. MODEL BOUNDARY LAYER TREATMENT 

The model BL parameters in regions of strong sensible 
heating are characterized by an oscillatory tendency which 
becomes apparent as the BL is iteratively developed at each 
evaluation time. COMP3 is successively iterated at each time 
to allow for boundary layer development, as previously dis- 
cussed. The BL parameters of thickness, and inversion jumps 
of wind, temperature and moisture are checked for convergence 
at each iteration. Oscillations, especially in the BL thick- 
ness, occur at points where stratus is diagnosed in the lowest 
layers. BL parameters at non-stratus gridpoints demonstrate 
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convergence within ten iterations. BL thickness oscillations 
with magnitudes in excess of 100 mb are observed at some grid- 
points. Similar oscillations have been previously noticed 
in the operational version of the model and a minor code 
change was instituted to reduce the effect. Discussions with 
Mr. Steve Payne of the Navy Environmental Prediction Re- 
search Facility (NEPRF) and further investigation provide 
the following scenario for model behavior at gridpoints where 
oscillations in the BL thickness fields are observed. 

The top of the boundary layer is restricted below the 
lowest sigma level (maximum BL thickness of about 200 mb) . 

A layer cloud instability (LCI) test is performed at grid- 
points where stratus is diagnosed by the model. Only stable 
stratus is permitted in the boundary layer. Unstable stratus 
occurs in areas where strong fluxes at the surface result in 
increased entrainment of dry air into the inversion region. 

This increased entrainment tends to eliminate the stratus deck, 
which is then treated as cumulus by the model. Since cumulus 
is not allowed within the boundary layer, the top of the 
boundary layer is lowered to the bottom of the cloud deck at 
gridpoints where unstable clouds are diagnosed. The moisture 
that was in the upper portion of the collapsing boundary layer 
is left behind to be removed as layer six large-scale 
precipitation. 

Oscillations in the BL thickness field are a result of 
the above process. As the boundary layer is developed 
through successive iterations at a gridpoint where unstable 
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stratus is identified, the BL height is allowed to increase 
until the unstable stratus is located within the boundary 
layer. The top of the boundary layer is then collapsed to 
the base of the unstable cloud- deck. In the next iteration, 
the boundary layer is allowed to grow again since unstable 
clouds are no longer within it. The following iterations will 
repeat this pattern and produce a bimodal oscillation in the 
BL thickness field. 

Similar boundary layer behavior was observed by Elsberry 
et al . (1984) in an evaluation of the NOGAPS diagnosis of the 

diabatic processes in eastern Pacific Ocean post-frontal 
convective clusters. That study emphasizes the modification 
of the transfer coefficients used in the model surface flux 
parameterization which occurs when the BL height decreases 
to a minimum value. The model sensible heating parameteriza- 
tion is given by Deardorff (1972) . His surface heat flux 
relationship is 



h = pc ( e - 0 ) u c C Q 

p s m m u 6 



( 2 ) 



where C is the friction transfer coefficient, C Q is the 
heat transfer coefficient, 0 g is the surface potential 
temperature, 6^ is the mean BL potential temperature, 
is the mean BL wind component, p is the density, and 
is the specific heat at constant pressure. The coefficients 
and C p are increased at gridpoints where the boundary 
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layer thickness becomes extremely small (< 250 m) . The very 
large surface sensible heat flux maxima produced by NOGAPS 
may then be due, in part, to the increase of the transfer 
coefficients caused by the collapsing boundary layers in 
these areas. 

The code change implemented in the operational version 
of the model, provided by Mr. Steve Payne of NEPRF, consists 
of an artificial boundary layer inversion which is inserted 
when LCI is diagnosed. This change reduces the magnitude of 
the oscillations that occur. A similar result was documented 
by Elsberry et al. (1984) when this change was implemented 
in their diagnosis. 

All NOGAPS fields are averaged in time over the last six 
iterations to give more representative results. This 
averaging is necessary because the implemented code modifica- 
tion leads to some trimodal as well as bimodal variations 
of reduced magnitude. The first four iterations are excluded 
from the averaging process to allow for initial convergence 
of the BL parameters at stable gridpoints. 

To verify model behavior in regions where instabilities 
are expected, fields of sensible heat flux, LCI, BL thickness 
and layer six large-scale precipitation are compared at each 
evaluation time. Fields from 12 GMT 13 January 1979 are 
shown in Figs. 2 and 3 for illustration purposes and are 
representative of other time periods. Figs. 2a and 2b present 
the surface sensible heat flux field and areas where LCI is 
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diagnosed by the model. A comparison at locations A and B 
of these figures clearly demonstrates that the model diagnoses 
unstable stratus and strong positive sensible heat flux in 
the same region. 

Locations of the small (20 mb or less) BL thickness and 
layer six large-scale precipitation are mapped in Fig. 3. 

A comparison at locations A and B of Figs. 3a and 3b with 
those of Figs. 2a and 2b shows that areas of minimum BL 
thickness and layer six large-scale precipitation are asso- 
ciated with the sensible heat flux maxima and LCI fields. 
Regions of large-scale precipitation not associated with the 
layer cloud instability are assumed to be the result of large- 
scale processes only. 

Comparisons of these figures support the scenario of 
model behavior as discussed previously. Evidence indicates 
that the model diagnoses unstable clouds where the sensible 
heating maxima occur. In these regions, the BL heights have 
collapsed, and moisture is removed as large-scale precipita- 
tion. This behavior illustrates an elaborate method for the 
model to account for the presence and effect of cloud streets 
and open-cell cumulus which form in the cold air where the 
surface sensible heating is strong. Although open-cell 
cumulus is actually a convective cloud, it is not labelled 
as such in the model. 
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D. EVALUATION OF SENSIBLE HEAT FLUX AND CLOUD FIELDS 

The NOGAPS fields of surface sensible heat flux, convec- 
tive and large-scale cloud location, and convective and 
large-scale condensation heating are evaluated through 
comparisons with ship reports, ECMWF FGGE analysis, satellite 
imagery and past studies. The purpose of the evaluation is 
to determine the validity of the fields as produced by the 
model. The 1000 mb wind analyses are produced directly from 
the FGGE wind components, while the diabatic heating, boundary 
layer thickness and layer cloud instability fields are taken 
directly from the model output. Convective and large-scale 
cloud maps are generated by assigning a value of one to each 
gridpoint where that type of precipitation is present in the 
model output, and a zero otherwise. The values at the six 
sigma levels are then summed at each time and the resulting, 
fields are plotted. Contours of these composite fields indi- 
cate where, and at how many levels, convective or large-scale 
precipitation is present in the model. 

1 . Sensible Heat Flux 

The evaluation of the sensible heat flux is based on 
the bulk formulation where sensible heating is directly 
proportional to the air-sea temperature difference and the 
boundary layer wind speed at each gridpoint. The 1000 mb 
ECMWF wind field is taken as representative of the boundary 
layer wind field for this evaluation. The NOGAPS surface 
temperature and BL surface air temperature, which is a model 
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extrapolation of the 1000 mb FGGE air temperature, are used 
to compute the air-sea temperature differences. Figures of 
satellite imagery are referred to out of sequence in this 
section to support the sensible heat flux diagnoses without 
complicating the comparison with the model cloud maps in the 
next section. 

a. 00 GMT 13 January 1979 

The storm system organized during the previous 
24 h, but has not reached the explosively deepening stage. 

The center of the disturbance is located approximately 600 
n mi to the southeast of Japan on the warm side of a relatively 
weak sea-surface temperature gradient. An intense cyclone is 
located to the east of the Kamchatka Peninsula. This cyclone 
dominates the flow pattern in the boundary layer wind field 
as illustrated in Fig. 4a. 

Two distinct features are evident in the sensible 
heat flux field (Fig. 4b) . The first is a two-pronged posi- 
tive flux maximum in the vicinity of, and to the north of, 
the storm center. A comparison with Fig. 4a indicates that 
this flux feature is tied closely to the BL wind field. The 
northern area, identified as 'A', results from the strong 
westerly flow to the southwest of the Kamchatka cyclone. 

This flow carries cold, continental air over a relatively 
warm sea surface, although the temperature gradient is weak. 

The magnitude and extent of this flux maximum is directly 
linked to the strength of the flow in the region, which peaks 
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at 27 m/sec, and the air-sea temperature difference, which is 
strongly negative over the area (-11 to -23°C) . Verifying 
ship reports from this region generally support the range 
of air-sea temperature differences determined from the ECMWF 
analyses (observed differences: -9 to -20 °C) . However, the 

observed windspeeds are of only half the magnitude of those 
represented in the analyses. 

The southern maximum area, identified as feature 
'B' in the figure, results from a separation in the westerly 
flow to the north, which transports some of the cold flow 
southward across the sea-surface temperature gradient. This 
split in the flow pattern is a manifestation of the develop- 
ing circulation of the incipient cyclone to the south. The 
wind speeds are generally lower, and the air-sea temperature 
difference (maximum: -17°C) is less than in the area of the 

northern maximum, which accounts for the smaller sensible heat 
flux diagnoses in this region. In situ measurements support 
the air-sea temperature difference and the windspeeds used 
in the model diagnoses for this region. 

The second feature of interest in Fig. 4b is the 
maximum of negative flux, identified as 'C', found to the 
east of the developing storm center. Although sizeable nega- 
tive surface fluxes are rare in a mid-latitude winter situa- 
tion, this negative flux is driven in the model by the strong 
southerly flow which funnels warm, tropical air across the 
sea-surface temperature gradient over colder water. A 
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comparison of the figures indicates that the region of nega- 
tive flux is coincident with this southerly flow. The 
largest negative flux is found where the magnitude of the 
southerly flow is highest (27 m/sec), and the air-sea 
temperature difference is most positive (+15°C) . Ship measure- 
ments from this region support an air-sea temperature differ- 
ence of only +2°C, and a maximum windspeed of 15 m/s. It 
appears that discrepancies in the ECMWF analyses force the 
model to diagnose excessive amounts of negative sensible heat 
flux . 

A comparison of Fig. 4b with the visual imagery 
for this time (Fig. 9) indicates that the region of sensible 
heat flux maximum to the north (A) has an extensive cover of 
cloud streets and open-cell cumulus which verifies the 
presence of strong surface fluxes, while sensible heat flux 
maximum 'B' is covered by the main storm cloudiness. Feature 
' C' appears to lie to the east of a zone of frontal cloudiness 
in the imagery, which would place it in the warm sector where 
low-level stratus would be expected, 
b. 12 GMT 13 January 1979 

The storm center moves northeastward over the 
largest sea-surface temperature gradient and enters the explo- 
sive stage. The general flux pattern (Fig. 5b) is similar 
to that of 00 GMT 13 January 1979. Maximum features in the 
positive flux field, identified as 'A' and 'B' in the figure, 
have become more distinct. This change in the flux pattern 
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occurs as the circulation associated with the developing 
cyclone intensifies and begins to dominate the flow pattern 
to the southwest, as depicted in the 1000 mb wind field 
(Fig. 5a). Verifying ship reports are not available for 
this time. 

Positive flux feature 'A' has elongated from 

west to east as the westerly flow to the south of the 

Kamchatka low now extends farther to the east. This cold, 

continental flow over the ocean maintains at least a -5°C 

air-sea temperature difference to 180 °E. Again, the highest 

flux values are found where the combination of high wind 

speed and large negative air-sea temperature differences 

are located. The magnitude of the maximum which stretches 

2 

to nearly 170 °E (600 W/m ) may be an overestimate through 

alteration of the drag coefficients by the oscillating 

boundary layer in this area. 

A comparison with the sensible heat flux field 

from the previous time indicates that the model correctly 

diagnoses smaller fluxes where this westerly flow weakens 

along the coast of Sahkalin Island. This coastal maximum 

2 2 

is decreased from approximately 1300 W/m to 900 W/m due 

to a 10 m/sec to 20 m/sec slackening of the wind in the region, 

while the air-sea temperature difference remains unchanged. 

Feature 'B' of Fig. 5b has also become elongated 
in the east-west direction and the maximum flux value has 
decreased in magnitude. The northerly flow, which breaks 
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off from the cold airstream to the south of Kamchatka, 
broadens to the east as more of this flow is pulled into the 
vicinity of the storm center by the increased storm circu- 
lation illustrated in Fig. 5a. This northerly flow brings 
cold air due south across the sea-surface temperature gradi- 
ent. The effect of this flow on the magnitude of the flux 
feature is diminished by the intrusion of tropical air to the 

east of the storm center. The location of the flux maximum 
2 

of 500 W/m to the east of Japan is coincident with the posi- 
tion of the maximum air-sea temperature difference in the 
area (-15°C) . The southern finger of feature 'B' is a mani- 
festation of the colder, westerly flow which passes to the 
south of Japan with a windspeed maximum in excess of 27 m/sec. 

The positioning of the negative flux feature ' C' 
in the figure is consistnet with the previous period. Strong 
southerly flow is still indicated in Fig. 5a to the east of 
the developing cyclone, but the magnitude and extent of the 
windspeed maximum has decreased and its position has shifted 
to the northeast. The negative sensible heat flux maximum 
is decreased from the previous time and its position is 
shifted to the northeast, which is consistent with the 
changes in location and magnitude of the wind speed and 
positive air-sea temperature difference. 

Visual satellite imagery from 1405 GMT 13 January 
(Fig. 12) clearly indicates the presence of cloud streets 
and open-cell cumulus in the cold air streams over both 
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positive sensible heat flux features in Fig. 5b. The nega- 
tive flux feature to the southeast of the storm center is 
properly located in the warm sector to the east of the main 
frontal cloudiness depicted in the imagery, 
c. 00 GMT 14 January 1979 

The storm center continues to move northeastward, 
but is still located over the region of strongest sea-surface 
temperature gradient. Explosive cyclogenesis continues and 
the cyclonic circulation increases, with a wind speed maximum 
of 30 m/sec (Fig. 6a) . 

The diagnosed flux field (Fig. 6b) merges the two 

major positive features (A and B) as the cyclone moves toward 

the dominant Kamchatka system and its circulation joins in 

the cold flow to the west of the developing storm center 

2 

(Fig. 6a) . The 800 W/m maximum in feature 'B ' (in the 
vicinity of 45°N, 165°E) is driven by the strong (15 m/sec) 
northerly flow of the storm circulation. The presence of 
this flow in an area with an air-sea temperature difference 
of approximately -15 °C supports the location of the flux 
maximum, although the magnitude is probably too high as noted 
previously. Ship reports from this region indicate that the 
observed air-sea temperature difference was approximately 
-12 °C and the wind speed was 12 m/s, which supports the 
ECMWF analyses. 

Feature 'A' now extends across the North Pacific 
to Alaska, and accompanies the westerly flow which is present 
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in this region. The magnitude of this flux is decreased to 
the east, which reflects the diminished air-sea temperature 
difference as the cold air has been modified by the long 
path length over the ocean. It seems that the magnitude of 
this flux feature has been overestimated to the east of 180 °E, 
but the moderate winds (10 m/sec) and relatively large air- 
sea temperature difference (-8°C) which extend into this 
area support the model flux diagnosis. Ship reports in the 
area generally support the computed air-sea temperature 
difference (observed: -4° to -8.3°C) and ECMWF wind speed 

(observed: 6 to 15 m/s) . 

2 

A negative maximum of approximately 300 W/m , 
identified as feature 'C' of Fig. 6b, is located on the eastern 
side of the developing cyclone. As in the previous cases, 
a strong southerly flow is driven by the large-scale circu- 
lation. This flow is located between 160°E and 170°W in the 
figure. It carries warm, moist tropical air over the colder 
ocean in the vicinity of the storm center, and results in 
positive air-sea temperature differences. Comparison of the 
figures shows that the negative flux maximum is located in 
the same position as the highest wind speeds (30 m/sec) over 
an area of positive air-sea temperature difference. It 
appears that the diagnoses of the presence of this feature 
is consistent with the analyzed low-level circulation for 
this time. No ship data are available for comparison in 
this area. 
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Visual satellite imagery from 0023 GMT 14 January 
1979 (Fig. 15) shows extensive open-cell cumulus coverage 
above 4 5°N. This region corresponds to the upper portion of 
the positive flux feature in Fig. 6b. More open-cell cumulus 
is indicated in the imagery to the west of 165°E between 
30 °N and 35°N, which corresponds to the southern portion of 
the positive flux feature. The flow pattern in the vicinity 
of sensible flux maximum 'B' (40°N, 160°E) is distorted by 

a polar low which has propagated into the area from the west. 
However, it appears that open-cell cumulus verifies in the 
imagery at this location as well. IR imagery (Fig. 16) 
indicates the presence of closed-cell cumulus in the southern 
portion of negative flux feature 1 C' to the east of the main 
frontal clouds which is consistent with the negative flux 
diagnoses. However, it appears that most of the region of 
diagnosed negative flux is covered by the main storm 
cloudiness . 

d. 12 GMT 14 January 1979 

The cyclone continues to grow explosively, and 
reaches peak intensity near this time. Fig. 7a illustrates 
a markedly different 1000 mb wind pattern to the north of the 
cyclone center from that previously analyzed. The westerly 
cold stream which has persisted to the south of Kamchatka 
is interrupted by the northward movement of this well- 
developed-cyclone, which is now positioned on the cold side 
of the sea-surface temperature gradient. Consistent with 
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the absence of the westerly cold stream to the north of the 
storm center at 1000 mb, the model diagnoses minimal sensi- 
ble heat flux in this region, as illustrated in Fig. 7b. 

A weak, positive flux feature, identified as 'A' in Fig. 7b, 
is located to the west of the storm center. This maximum 
results from the 15 m/sec to 25 m/sec winds which separate 

from the relatively cold, southwesterly flow at about 48°N, 

2 

170°E. The positive 500 W/m maximum to the south of Kamchatka 
is driven by the continental flow as it departs the coast at 
this higher latitude. 

Feature 'B' of the figure identifies a very large 

2 

and unrealistic negative flux maximum of more than 700 W/m 
to the southeast of the storm center. This anomalous feature 
appears to originate with the ECMWF analyses rather than model 
error. Excessive windspeeds of greater than 50 m/sec have 
been analyzed in the southerly flow in this region. The 
model diagnosis seems consistent with this large magnitude of 
tropical flow in an area of positive air-sea temperature 
difference (8°C) . Ship data are not available from this 
area for comparison. 

Open-cell cumulus verifies in the visual imagery 
from 1205 GMT 14 January 1979 (Fig. 18) for the region of 
positive flux feature 'A' in Fig. 7b (65°N-45°N, 160°E-175°E) . 
The IR imagery (Fig. 19) shows extensive closed-cell cumulus 
and stratus coverage in the warm air to the east of the frontal 
zone above 20 °N. This area generally corresponds to the 
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region of negative flux diagnosed by the model. However, 
most of this diagnosed negative flux feature is covered by 
the main storm cloudiness. 

e. 00 GMT 15 January 1979 

The storm has occluded over the cold water in 
the vicinity of the Aleutian Islands at this time. A cyclone 
which has formed to the southwest of the storm center, and 
an anticyclone which is located to the southeast, tend to 
force a southerly flow into the storm region, as illustrated 
in Fig. 8a. 

2 

A negative flux feature of approximately 200 W/m 
is located to the southeast of the storm center. This feature 
is a consequence of the strong southwesterly * flow (35 m/sec) 
to the southeast of the storm center, which is forced by the 
circulation of the incipient system and the large anticyclone 
to the southeast. Verifying ship observations are not avail- 
able from this area. A comparison of the figures indicates 
that this maximum is located where the southerly flow splits 
just to the southeast of the storm center, in the region 
where maximum wind speeds of 20 m/sec to 25 m/sec occur. The 
positive air-sea temperature difference has been decreased 
due to the modification of the air as it reaches this high 
latitude. The presence and extent of this negative flux 
feature appears to be dominated in the model by the magnitude 
of the wind. 

The location of the positive sensible heat flux 
feature to the west and south of the storm center, identified 
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as 'A' in Fig. 8b, is determined by the southwesterly and 
westerly flow in the vicinity of 50°N. This moderate, modified- 
continental flow serves as the source of the cold air to the 
west of the storm center and maintains a small negative air- 
sea temperature difference in the region. The small tempera- 
ture difference and the general weakness of the flow leads 

2 

to a model diagnosis of fluxes of 100 W/m or less to the 
west of the storm center. It is interesting to note that 
this feature begins to wrap around the storm center to the 
south for the first time. The intrusion of this feature into 
this area of the storm is indicative of the occlusion process 
which is taking place at this time. 

f. Sensible Heat Flux Summary 

The model boundary layer parameterization leads 
to diagnoses of sensible heating that are consistent with 
the forcing represented in the ECMWF analyses. The positive 
fluxes are confined to the cold air behind the front, which 
conforms to Petterssen's general cyclone model, and negative 
fluxes are located in the warm, tropical air to the east. 

The flux maxima are consistently placed where the greatest 
air-sea temperature difference and boundary layer windspeed 
coincide . 

The magnitudes of the fluxes generated by the 
model seem to be excessive. The amount of positive flux is 
affected by discrepancies in the ECMWF analyses, as well as 
the artificial increase in the drag coefficients where the 
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boundary layer is collapsed by the model. The unrealistic 
magnitudes of the negative fluxes present in some of the 
fields appear to be related to spurious ECMWF 1000 mb wind 
analyses and the large positive air-sea temperature differ- 
ences determined in these areas. 

Maximum sensible heat flux values found by other 
investigators are included in Table I for comparison purposes. 
These values were compiled by Gall and Johnson (1971) with 

an addition from Bosart and Lin (1984). Units have been 

2 

converted to W/m where applicable. The magnitude of the 
flux produced by NOG.APS is of the same order as that deter- 
mined by previous investigators, and may be representative 
of the intense nature of this cyclone. 

2 . Cloud Verification 

Storm-related composite cloud features, as generated 
by the model, are compared to Defense Meteorological Satellite 
Program (DMSP) imagery at each analysis time. The discussion 
focuses on three general model cloud forcing categories. 

These categories are convective, large-scale and cumulus, 
where cumulus identifies the open-cell cumulus and cloud 
streets that are diagnosed as layer-six, large-scale clouds 
in the model. Cumulus occurs where the sensible heat flux 
maxima are located, as discussed in Section C of this chapter, 
a. 00 GMT 13 January 1979 

Visual and infrared (IR) imagery for this time 
are included as Figs. 9 and 10. The composite convective 
and large-scale clouds are mapped in Fig. 11. 
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(1) Convective Clouds. The major cloud feature 



associated with the system at this time is identified as 'A' 
in the figures. The convective cloud map (Fig. 11a) denotes 
a narrow band of convection which penetrates three sigma levels 
to the west of the storm center. A comparison of the visual 
and IR satellite imagery shows that significant convection 
did occur in the vicinity of the storm center at this time. 
However, much of this convection is located to the east of 
the center in an area roughly bounded by 32 °N to 41°N and 
150 °E to 160 °E. The clouds that are present to the northeast 
of the storm center in the imagery are not diagnosed as 
convective by the model . 

An extensive band of convection is indicated 
to the east of the storm and is identified as feature 'B' 
in the figures. The IR imagery of Fig. 10 indicates that 
deep convection with low cloud-top temperatures occurs along 
the entire length of this band. The model correctly diag- 
noses deep convection which extends through four sigma 
levels as illustrated in the cloud map, but displaces it 
approximately 10 degrees to the west of the feature present 
in the imagery. The displacement of this feature is consis- 
tent with that observed in feature 'A'. 

(2) Large-Scale Clouds . The IR imagery indi- 
cates that low-level clouds are present to the east of the 
storm center in the same region where large-scale cloudiness 
is diagnosed by the model (centered at approximately 34 °N, 
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156°E) . The three- level, large-scale cloud feature located 
north and northwest of the storm center appears to be repre- 
sentative of the clouds present in this region of the satel- 
lite imagery. This cloud location with respect to the storm 
center infers large-scale cloud development as the tropical 
flow, now involved in the developing cyclone circulation, 
overrides the colder air to the north. The depth of the clouds 
generated by the model in this area is consistent with the 
moderately low cloud-top temperatures indicated in the IR 
imagery. The model diagnoses of the large-scale process in 
the warm air to the east and northeast of the storm center 
appear to be representative, although the area coverage of 
the large-scale clouds observed in the imagery is more 
extensive than produced in the model. This difference in 
coverage may be due, in part, to inactive clouds which are 
present in the imagery, but can not be diagnosed by the 
model . 

(3) Cumulus . The westward and northward extent 
of feature 'A' of Fig. lib can not be verified by the imagery, 
but occurs where the cold air separates from the westerly 
flow to the north, and is caught in the circulation of this 
system. This portion of the large-scale cloud feature is 
co-located with an area of sensible heat flux maximum, and is 
more than likely representative of the open-cell cumulus. 

The extensive region of single level large-scale clouds to 
the north of the storm center (above 40°N) in the figure is 
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generated under similar conditions. In this area, the visual 
imagery clearly shows the presence of the open-cell cumulus 
in the cold, westerly flow streaming off Sahkalin Island. The 
presence of this cumulus in an area where low-level large- 
scale clouds are indicated supports the scenario of the 
boundary layer behavior as described in Section C. 
b. 12 GMT 13 January 1979 

The visual and IR satellite imagery for this time 
are shown in Figs. 12 and 13. Convective and large-scale 
composite cloud maps are displayed in Fig. 14. A more 
extensive storm-related cloud mass is present in the imagery 
as the cyclone enters the explosive stage. 

(1) Convective Clouds . The satellite imagery 
depicts an extensive, comma-shaped cloud pattern associated 
with this cyclone, with significant convection occurring in 
the vicinity of the storm center. The tail of the comma 
extends from about 35 °N southeast to 20 °N. The IR imagery 
indicates that deep convection occurs in a narrow band, 
roughly 120 n mi in width, along the front edge of the tail. 
The convective cloudiness, as mapped by the model, displays 
a relatively deep and extensive pattern in the vicinity of 
the storm center. Portions of this model feature to east 
of 165°E extend beyond the coverage of the imagery. 

A comparison of the imagery and the convec- 
tive cloud map, Fig. 14a, reveals that the moderate storm 
convection to the west of the center, identified as feature 
'A' in the figures, is represented accurately by the model. 
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As in the previous time, convection is not indicated to the 
northeast of the center where significant deep cloud cover- 
age is shown in the imagery. The model also overestimates 
the extent of' the convective activity in the frontal band. 

The model diagnoses convection (Fig. 14a) extending through 
four sigma levels over a much broader region than can be 
supported by the imagery in this zone. The positions of 
areas of deep convection that are displayed in the imagery 
and cloud map correspond well. This comparison indicates 
that the model diagnosis is basically correct, although the 
modeled clouds are extended too far to the east. The exten- 
sive area of convective coverage farther to the south is 
considered a model anomaly. Overestimation of convection 
by the model along the southern border is present at each 
of the subsequent times. This behavior is described in more 
detail in Subsection D.3 of this chapter. 

(2) Large-Scale Clouds . Model-diagnosed, large- 
scale clouds are present in one to three sigma levels to the 
east and northeast of the storm center. This position and 
area of coverage are consistent with what is observed in the 
imagery for this region of the cyclone. Low cloud-top tem- 
peratures inferred from the IR imagery indicate that these 
clouds are not shallow, which supports the multi-level model 
diagnoses. As noted previously, the model correctly asso- 
ciates this cloudiness with the large-scale tropical flow in 
the region. When this cloud feature is combined with the 
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convective map, a representative depiction of the cloudiness 
as observed in the imagery is produced. 

(3) Cumulus. The single-level, large-scale 
cloudiness to the west and northwest of the storm center 
(Fig. 14b) corresponds to the location of the sensible heat 
flux maximum in the storm area. The visual imagery shows that 
open-cell cumulus is present in the cold air streaming behind 
the front in this area. A similar open-cell cumulus pattern 
is indicated to the north of this system where single-level, 
large-scale clouds are mapped over the sensible heat flux 
maximum in the westerly flow off Siberia. The diagnosis in 
these areas where the oscillating boundary layer behavior was 
observed is consistent with the expected treatment of the 
open-cell cumulus by the model. 

c. 00 GMT 14 January 1979 

The visual and IR satellite imagery is included 
as Figs. 15 and 16. The convective and large-scale composite 
cloud maps are included as Fig. 17. The satellite imagery for 
this time illustrates the continued development of this sys- 
tem, with enhanced cloud cover extending above 40°N. 

(1) Convective Clouds . As in the previous 
analysis times, the region of most extensive deep cloud 
coverage to the northeast of the storm center is not iden- 
tified as convective. The cloud pattern generally maintains 
a comma shape with more organized deep convection to the west 
and northwest of the storm center (Fig. 16) . A much broader 
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band of deep convection is present in the frontal zone ex- 
tending to the southwest. The structure of this frontal con- 
vection is diagnosed correctly by the model. A comparison 
of Fig. 17a with the IR imagery shows that the model provides 
a representative display of the orientation, depth and extent 
of the frontal convection as it occurred above 30 °N. However, 
the overall feature is displaced between five and ten degrees 
to the west of the frontal feature observed in the satellite 
imagery. A similar displacement is also noted in the 00 GMT 
13 January 1979 convective field. It appears the convective 
process is overestimated to the south of 30 °N. 

The small disturbance which is entering the 
storm region from the west is identified as feature 'A' in 
the satellite imagery. This disturbance may be captured in 
•feature 'A' of the model diagnosis (Fig. 17a). The three- 
level convection indicated by the model is consistent with 
the moderately low cloud-top temperatures inferred from the 
IR imagery. The leading edge and northward extent of this 
feature are well positioned by the model, but the overall 
area coverage is too extensive. A second possibility is 
that this diagnosed convective cloud area is a western dis- 
placement of convection associated with the major cloud system. 

(2) Large-Scale Clouds . The diagnosis of the 
large-scale clouds (Fig. 17b) is similar to the previous 
analysis time. Multi-layered large-scale clouds are indicated 
to the north and northeast of the storm center in the 
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satellite images. The model diagnoses large-scale clouds 
at up to four sigma levels where some of the deepest clouds 
are present in the satellite imagery. The single-level, 
large-scale feature, identified as 'B' in the figures, 
corresponds to the position of the low-level stratus present 
behind the frontal band (see especially Fig. 16). 

(3) Cumulus . The single-level, large-scale 
cloud feature to the west and northwest of the storm center 
in Fig. 17b is coincident with the sensible heat flux maxi- 
mum diagnosed at this time. The visual imagery indicates the 
presence of some open-cell cumulus to the west of the storm 
center and an abundance of this cloud type to the northwest. 
This cumulus distribution falls within a region diagnosed by 
the model as having low-level, large-scale clouds. This 
treatment of the open-cell cumulus is consistent with the 
previous evaluation times. 

d. 12 GMT 14 January 1979 

The visual and IR satellite imagery for this 
evaluation time are given in Figs. 18 and 19. The convective 
and large-scale composite cloud features are mapped in Fig. 
20. The visual imagery indicates that the cyclone is well- 
developed. The storm center is now completely surrounded 
by an extensive cloud mass except for an area of cold air 
penetration from the south. 

(1) Convective Clouds . The IR satellite imagery 
displays moderate convection in the vicinity of the storm 
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center. Deep convection occurs to the east and along two 
bands in the frontal zone, where it appears that a secondary 
front has formed. The disturbance that was previously noted 
to the west has developed considerably during the previous 
12 h. This disturbance is now located to the southwest of 
the storm center and is indent if ied as feature 'A* in the 
figures. A comparison of the imagery and the convective cloud 
map again shows that the deep clouds to the east and north- 
east of the storm center are not identified as convective 
by the model. The convective feature that is associated with 
this cyclone and that of disturbance 'A' have been merged 
by the model to the west of the storm center (Fig. 20a) . 

The east-west and north-south extent of feature 'A' has 
been exaggerated in the model in comparison with the satellite 
imagery. The position of the frontal band (Fig. 20a) corres- 
ponds very closely to that of the secondary front in the IR 
imagery. The band of deep convection to the east of this 
frontal feature is not represented in the model display. 

(2) Large-Scale Clouds . The model-generated, 
large-scale cloud map again complements the convective map. 

As in previous times, multi-layered, large-scale clouds are 
positioned to the northeast of the storm center. However, the 
diagnosed area coverage of the cloudiness in this region is 
much less than observed in the satellite imagery. The location 
of the large-scale clouds in the frontal band (Fig. 20b) pro- 
duces a very representative display of the position of the 
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observed cloudiness, especially in the eastward extent of 
the frontal band. For the first time, multi-layered (three 
level) large-scale clouds are located in the frontal zone. 

The presence of these clouds in the model is supported by 
the deep clouds that are located in this region of the IR 
imagery, and may be indicative of greater large-scale frontal 
activity as the storm reaches its peak intensity. Most of 
the frontal clouds are diagnosed as single-level however, 
and can not account for the deep penetration indicated in the 
IR imagery to the east of the storm center. Low-level clouds 
are present throughout the frontal band in the IR imagery, 
as well as in the warm sector to the east. The broadening 
• of the large-scale frontal feature in the model cloud map 
to the north of 40 °N appears to capture some of this warm 
sector stratus. 

(3) Cumulus . The single-level, large-scale cloud 
feature, identified as 'B 1 in Fig. 20b, corresponds to the 
position of the positive sensible heat flux maximum. The 
location of the open-cell cumulus, as shown in the imagery, 
is primarily to the northwest of 50°N, 170°E. Only the 
northwest corner of feature 'B' verifies as open-cell cumulus, 
with the remainder of the feature occupied by the main cloud 
cover of the cyclone. The open-cell cumulus to the north of 
55°N in the imagery also is not captured in the model diag- 
noses. The small finger of feature ' B' which extends to the 
east and connects the body of the feature with the storm 
center, is not located in an area of positive flux. The IR 
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imagery indicates that the clouds in the cold air in this 
area are limited to the low level and are probably large- 
scale in nature. Convection would be impeded in this area 
by the presence of the positive air-sea temperature 
difference . 

e. 00 GMT 15 January 1979 

Satellite imagery for this time (not shown) , 
and the previous analysis of this system by Calland (1983) 
shows that the cyclone has entered the occluded stage. 

(1) Convective Clouds . The convective cloud map 
(Fig. 21a) indicates that no storm-related convective clouds 
are diagnosed by the model. The large area of convection 
which is mapped to the west of the storm center is associated 
with the second disturbance previously mentioned. The diag- 
nosis of no convection at this time seems extreme in light of 
satellite imagery from this period and later periods which 
indicates that the system retains a considerable cloud pat- 
tern. However, this diagnosis is consistent with the occluded 
nature of the cyclone. The storm has moved over much colder 
water and the convective process has been effectively shut 
off in the model. 

(2) Large-Scale Clouds . The large-scale cloud 
features (Fig. 21b) do indicate that extensive cloud cover 
is still associated with this system. This large-scale pat- 
tern about the storm center includes multi-layered clouds 

to the north and northwest, and a blurred, single-level 



frontal zone curving to the southwest. This pattern and 
general coverage are very similar to that observed in the 
satellite imagery for 1005 GMT 15 January 1979 (not shown) . 

(3) Cumulus . The model behavior in the region 
of positive sensible heat flux appears to be similar to that 
of the previous times. The large-scale cloud feature in 
the cold air to the south, west and northwest of the storm 
center is positioned where the maximum positive sensible heat 
flux has been diagnosed by the model. Open-cell cumulus 
verifies in this region of the imagery, 
f. Cloud Comparison Summary 

The model-generated convective cloudiness is 
generally representative of the storm environment as verified 
by the satellite imagery. The outbreak of convection in the 
vicinity of the storm center and the deep convection over 
the frontal surfaces are correctly diagnosed. However, 
there is a tendency to displace or extend these features 
several degrees to the west of the verifying position. 

The open-cell cumulus and cloud streets located 
in the cold air streams are treated as large-scale convection 
in the NOGAPS boundary layer parameterization. If this 
interpretation is accepted, the generation of this cloud 
type is handled well by the model . 

Large-scale clouds are consistently placed to the 
north and northeast of the storm center, as well as along the 
frontal surfaces. The most extensive cloud mass is observed 
in the imagery in these areas during the period of cyclone 
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growth. The identification of large-scale cloud forcing in 
these cyclone regions appears correct. When the model output 
from the three forcing categories is combined, a reasonably 
accurate depiction of the storm clouds which are present in 
the imagery is produced. 

3 . Cumulus and Large-Scale Cloud Model Heating Rates 

The rate of condensation heating in each model sigma 
layer is a sum of the heating from the existing cloud types. 
Precipitation rates are combined in the model to produce a 
total precipitation field for each evaluation time (not 
shown) . 

Temperature tendencies from latent heat release for 
the six sigma layers are averaged to produce a column tem- 
perature change from cumulus and large-scale precipitation. 

The layer temperature changes are weighted by the thickness 
of the individual layers prior to column averaging. It 
should be noted that the large-scale clouds are confined 
primarily to the lowest two sigma layers, while the deep 
convective clouds extend well into the upper layers of the 
atmosphere. Column averaging is more meteorologically signi- 
ficant for convective condensation heating, because of the 
greater vertical extent of this cloud type. Averaging of the 
large-scale heating field artificially spreads its low-level 
contribution over the column. However, column averages allow 
a consistent comparison between the two precipitation sources 
and indicates the relative contribution of large-scale conden- 
sation heating to the total column temperature change. 
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A general feature which is noticed at all times is 
the excessive convective heating rates that are generated 
along the southern border of the figures. Here the model 
apparently responds to the tropical air and produces large 
areas of anomalous deep convection. 

In his semi-prognostic test of the Arakawa-Schubert 
cumulus parameterization, Lord (1982) found that the model- 
generated cumulus warming of the column and associated 
precipitation rates closely matched observations. His pro- 
files of cumulus warming and time series of precipitation 
from cases of moderate tropical convection generally gave 
column-averaged layer heating rates of less than 5 deg/day 
and time-averaged precipitation rates of 10 mm/day to 20 
mm/day. Low-tropospheric values of cumulus warming ranged 
between approximately 3 deg/day to 6 deg/day. The maximum 
precipitation rates observed were about 35 mm/day. 
a. 00 GMT 13 January 1979 

The convective and large-scale condensation 
heating fields are illustrated in Fig. 22. At this early 
stage of cyclone development, the largest extent of heating is 
provided by the large-scale processes as shown in Fig. 22a. 

The low-level heating rates for this area are about 2.5 deg/ 
day, which compares favorably with the GATE observations used 
in Lord's study. 

Storm- area convective heating (Fig. 22b) is con- 
fined to a small area to the west of the storm center, with 
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a column-averaged rate of about 5-10 deg/day . This heating 
is also consistent with Lord's findings for cases of moderate 
convection. The large feature to the east of the storm 
center indicates convective heating in excess of 30 deg/day 
to the north of 30 °N. The convective heating rates to the 
south of this latitude are excessive and unrealistic, 
b. 12 GMT 13 January 1979 

A large-scale column average condensation heating 
maximum of 0.5 deg/day is indicated to the west of the storm 
center in Fig. 23a, which actually is associated with a 
maximum heating rate of 2.5 deg/day in the lowest layer only. 
The bulk of this heating is derived from the open-cell 
cumulus, which verifies for this area in the imagery. The 
3.0 deg/day heating feature to the east of the storm center 
reflects a single gridpoint diagnosis which has been contoured 
as an area. 

The storm center is adjacent to an area of convec- 
tive heating with a column-averaged maximum of 30 deg/day. 

This convective heating appears to be excessive, but is sup- 
ported by a surface precipitation rate of approximately 12 
cm/day. Although this rate is much higher than in Lord's 
study, it is not unbelievable given the explosive nature of 
the cyclone. An area of much greater heating is present in 
the feature below 28°N. This excessive heating is produced 
by the anomalous convection along the southern portion of 
the grid mentioned previously. 
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c. 00 GMT 14 January 1979 

Large-scale condensation heating occurs to the 
west, northwest and northeast of the storm center (Fig. 24a). 
The column-averaged heating rate maximum for the region to 
the west is approximately 2.5 deg/day. Since this heating 
area corresponds to the moderate convection of the open-cell 
cumulus in the cold air behind the storm center and to the 
north, it appears that the model calculations are of the 
correct magnitude. 

The convective precipitation heating rates (Fig. 
24b) are in excess of 30 deg/day in the frontal region of the 
cyclone. This heating rate seems excessive, as in the 
previous evaluation time, but it is supported by a diagnosed 
surface precipitation rate maximum of approximately 16 cm/d ay . 
This precipitation rate is almost 4.5 times greater than the 
maximum observed by Lord (1982). An increase in the heating 
rate can be identified in this field at about 31°N which 
separates the storm frontal convective heating from the region 
of anomalous heating to the south. 

d. 12 GMT 14 January 1979 

The maximum in the large-scale heating field 
(Fig. 25a) is located in the frontal zone to the southeast 
of the storm center. The actual low-level heating rate for 
this area ia 6 deg/day, which is reasonable given the inten- 
sity of the system. 

The magnitude of the convective heating is signi- 
ficantly reduced (Fig. 25b) over previous times. A distinction 
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can be made between the storm area heating and the excessive 
anomalous heating to the south in this field, as in the 
previous analyses. 

e. 00 GMT 15 January 1979 

The large-scale, column-averaged condensation 
heating rate is realistically diagnosed as 1.0-1. 5 deg/day 
in the frontal zone, and 0. 5-2.0 deg/day to the west and 
southwest of the storm center. The positioning of this 
large-scale heating properly accounts for the contribution 
of the open-cell cumulus in this area. 

No convective heating is indicated by the model 
in the vicinity of the storm. This behavior is a reflection 
of the occlusion process which has been correctly diagnosed 
by the model. The convective heating field (Fig. 25b) shows 
a region of intense heating to the south of the storm center. 
This heating feature can be attributed to the second distur- 
bance which has moved into the area, but may also include some 
of the anomalous heating that is observed at every other 
evaluation time. 

f. Condensation Heating Summary 

The column-averaged, convective condensation heat- 
ing rates in the frontal zones were found to be approximately 
30 deg/day, which appears to be high, but is not inconsistent 
with the surface precipitation generated by the model in the 
storm area. Surface precipitation on the order of 10 cm/day 
is not implausible given the intense nature of this cyclone. 
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Therefore, the large frontal convective heating rates diag- 
nosed by the model are probably representative of the FGGE 
forcing. Convective heating is definitely overestimated 
along the southern borders of the domain. It seems that the 
model is responding to both the tropical environment and the 
southern boundary in these regions. It is noted that NOGAPS 
does not see a boundary of this type in an operational run. 
The overestimation along the border therefore, is probably 
an artificial inducement not representative of NOGAPS 
performance . 

Large-scale heating is well represented by the 
model, especially where open-cell cumulus and cloud streets 
are present. In general, the low-level temperature change 
rates correspond to those reported by Lord (1982) , and the 
range of the column-averaged temperature change rates (0.5- 
2.5 deg/day) is physically reasonable for this category of 
cloud forcing. The largest contribution of large-scale con- 
densation heating was properly diagnosed in the frontal 
regions during the last two evaluation times when the cyclone 
reached peak intensity. Column- averaged heating rates in 
regions where cloud streets and open-cell cumulus verified 
in the imagery were consistent with Lord's results for cases 
of moderate convection. 



69 



IV. ADVECTIVE AND DIABATIC TEMPERATURE 
CHANGES DURING CYCLOGENESIS 



A. GENERAL 

The diabatic contribution to column temperature change 
in the storm environment is estimated through the application 
of a quasi-Lagrangian thermodynamic energy budget and from 
the NOGAPS diabatic diagnostics. The equation employed in 
the budget routine is derived from the first law of thermo- 
dynamics using the potential temperature as described by 
Holton ( 19 79 ) , 



_q_ 

c 

p 



,3 £n0 x Q ^ 

[ — + V • V £n 0 + a) 



3 £n 8 
3p 



] x T 



(3) 



In (3), q is the diabatic heating, C is the specific heat 

P 

at constant pressure, term 1 is the 12-hour time tendency 
of £n 0, term 2 represents the horizontal advection, term 3 
is the vertical advection, and T is the temperature in the 
layer. The units of q/C^ are K/day. 

The first procedure is to infer the right side of (3) 
using the quasi-Lagrangian budget technique as discussed by 
Calland (1983). The second procedure is to estimate the 
contribution of this heating term from the diagnostic re- 
sults of the NOGAPS heating package. 

The storm-area diabatic temperature changes are computed 
by layer and averaged to produce column diabatic temperature 
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change for the storm volume. Column-averaging permits a 
direct comparison of the budget results in isobaric coordi- 
nates with NOGAPS results in sigma coordinates, and yields a 
more accurate representation of the net effect of the thermal 
processes in the storm volume. Weighting factors are applied 
to layer heating values to normalize the results to a standard 
100 mb layer prior to column-averaging. The weighting fac- 
tors that are applied to the budget pressure layers and the 
NOGAPS sigma layers at 12 GMT 13 January are listed in Table 
II for illustration. Date-times that are included in the 
figures represent 12-h time periods. For example, 1206 
refers to 00-12 GMT 12 January 1979. 

B. THERMODYNAMIC ENERGY BUDGET CALCULATION 
1 . Procedure 

The quasi-Lagrangian thermodynamic energy budget 
program includes interpolation of the basic observed variables 
of temperature, horizontal velocity and vertical velocity 
to a given storm volume in increments of one degree radii 
from the storm center. Observed temperature tendencies, 
horizontal temperature advection and vertical temperature 
advection are determined from these interpolated fields. The 
budget terms are area-averaged by layer over the storm 
volume . 

a. Finite Differencing and Averaging 

Locations of the variables in the budget compu- 
tation are displayed in Fig. 27a. The budget terms are 
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structured to produce layer results. Temperatures and veloci- 
ties are averaged to the center of each layer except where 
vertical advection is determined. The vertical advection 
computation uses a layer-averaged omega, but carries the 
level temperatures to allow the finite difference (3 Jin 6/3P) 
to be centered in the layer. The final form of the budget 
equation is: 

— -A A A — A 

= (d IrTO/dt + V • V TnT + w 3 Jin 6/3P) * T (4) 

P 

where the terms are as defined previously and the overbars 
indicate layer and area averages, and V is the wind velocity 
relative to the storm. 

ECMWF temperature analyses are available at 00, 

06, 12 and 18 GMT for each day. Because of the possible 
errors in the 06 and 18 GMT potential temperature analyses 
(Calland, 1983) only the 00 and 12 GMT analyses are utilized 
in this study. Time resolution is sacrificed for added 
confidence in the temperature time tendencies using the 00 
and 12 GMT data. 

The time difference of Jin 0 is determined between 
adjacent observations as illustrated in Fig. 27b. The dia- 
batic temperature change results are determined from the 00 
and 12 GMT analyses and are valid for the 12 h period 
centered at 06 and 18 GMT due to the centered time differencing 
and the time-averaging of the advection terms. Vertical time 
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sections are constructed using the weighted results from 06 
and 18 GMT for the period 12-14 January 1979. 
b . Data 

Kinematic vertical velocities are used in the 
budget calculation because the ECMWF FGGE vertical velocities 
are not the product of direct analysis, but are determined 
during the initialization step of the 4 -dimensional data 
assimilation technique employed at ECMWF (Bengtsson, et al., 
1982) . The kinematic fields are chosen to ensure that the 
horizontal and vertical velocities, as well as the vertical 
finite differencing of the thermodynamic energy equation, are 
consistent with the continuity equation. The magnitude of 
the kinematically-derived omega is nearly two times that of 
the initialized FGGE omega. 

Mandatory-level temperatures are vertically 
averaged (linear on pressure) to obtain layer mean tempera- 
tures. These temperatures were determined during the 
initialization phase of the data assimilation cycle at ECMWF. 
To determine the magnitude of the differences between the 
initialized and analyzed thermal structures, the ECMWF layer 
temperatures are compared to equivalent layer temperatures 
computed by the hypsometric equation using the ECMWF height 
field. These heights are determined during the analysis cycle 
of the data assimilation. A comparison of these two tempera- 
ture fields lends confidence to the ECMWF level temperature 
data in the troposphere in the vicinity of the developing 
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cyclone. Temperatures determined by the two methods for the 
500-700 mb layer at budget radius 6 are displayed in Fig. 28. 
These temperatures are generally in agreement, but some 
significant ’ differences do exist. A 2.5 K difference is 
present in Fig. 28 at time period 6. Departures of this 
magnitude may adversely affect the budget calculation. 

The largest deviation between the two tempera- 
ture fields is found in the stratospheric layers, which sug- 
gests a problem with the ECMWF initialization in the upper 
three levels. The average deviation between the ECMWF and 
calculated temperatures in the troposphere is 0 . 84 K. The 
average stratospheric deviation is 2.73 K. The largest single 
departure in the troposphere is 3.0 K, while the maximum 
stratospheric departure is 7.3 K. The RMS difference between 
the initialized and calculated layer temperatures at radius 
6 for all levels at all times is 1.86 K. The initialized 
ECMWF temperatures were elected for the budget calculations 
to be consistent with the NOGAPS diagnosis, which was per- 
formed with this data. 

Since the analyzed temperatures are only carried 
at discrete levels, the model lapse rate will smooth out 
inversion features within the layers. This smoothed lapse 
rate could potentially introduce errors into the budget 
calculation, but the effect is greatly reduced through 
column- aver aging . 



74 



c. Heating and Temperature Change 

The budget residual represents a layer temperature 
change in units of K/day due to diabatic heating and accumu- 
lated errors in the budget calculation. This temperature 
change is expressed in per unit mass (Kg) and is proportional 
to the amount of heating (q) applied to the layer. These 
temperature changes cannot be compared directly because the 
layers in the column are of varying thickness (mass) . To 
compare temperature change rates from layers of unequal mass, 
all layer temperature change values are normalized to a 
standard layer of 100 mb thickness. The resulting weighted 
temperature tendencies are not indicative of the actual 
sigma or mandatory layer temperature changes which occur, 
but are proportional to the amount of heating that is applied 
vertically to the storm volume. 

Column-averaged temperature change rates are com- 
puted by vertically summing the weighted layer temperatures 
and dividing by the sum of the weighting factors. The result 
of this calculation is representative of the actual average 
column temperature change rate; assuming the budget calculation 
is free of errors. 

2 . Budget Analysis 

Vertical time sections of the results from the budget 
calculation are analyzed for the inner and outer storm 
volumes of six and ten degree radii, respectively. 
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a. Temperature Time Tendency 

The temperature tendency term for the inner volume 
(Fig. 29a) shows a temperature increase through the column to 
175 mb during the first two time periods. Most of the calcu- 
lated temperature increase is confined to the low levels with 
a maximum at 775 mb at 1218. A weaker area of upper-level 
increase is reflected during this period as well. This calcu- 
lated temperature increase precedes the onset of explosive 
cyclone development by approximately 12 h. The location and 
distribution of this temperature increase reflects forcing 
which will enhance cyclone intensification. 

Negative values dominate the tropospheric tempera- 
ture tendency during all periods after 1306. This feature is 
expected in the lower troposphere as the storm center moves 
northward into colder air. The zero contour which marks the 
top of this negative tendency feature can be treated as an 
indicator of tropopause behavior. The tropopause lowers 
between 1318 and 1406, indicating that the storm center 
transits into a colder air mass during these periods. The 
large negative maximum (-22.5 K/day) present at 1406 occurs 
concurrently with the movement of the storm across the 
atmospheric and oceanic baroclinic zone. This maximum re- 
flects the large low-level volume temperature change which 
occurs as the cyclone moves rapidly into the cold air. The 
location of this negative temperature tendency maximum at 
600 mb places it within the cold dome as well. 
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The vertical extent of the negative temperature 
tendency during the explosive period is inconsistent with 
continued cyclogenesis. The temperature changes represented 
in Fig. 29a suggest that the column contracted during the 
period of most intense development, when the opposite should 
have occurred. Although negative temperature tendencies were 
expected in the lower troposphere as the storm moved north, 
the lack of compensating warming in the upper levels is sus- 
pect. Previous evaluation of the initialized ECMWF temperature 
fields indicated that some discrepancies exist in the upper- 
level temperature analyses. These temperature discrepancies 
may have contaminated the budget tendency term in the upper 
levels . 

The temperature tendency time section for the 
outer budget volume (Fig. 29b) is similar to that of the 
inner volume. Many of the features have been smoothed and 
reduced in magnitude, which reflects the averaging over a 
larger region. 

b. Advective Temperature Tendency 

The advective tendency term at the inner volume 
(Fig. 30a) generally indicates a strong temperature decrease 
in the lower troposphere with a weak increase in the upper 
layers. Initial upper-layer positive temperature advection 
is centered at 275 mb. Horizontal warm advection dominates 
the upper-layer advective temperature tendency until explo- 
sive development becomes established after 1306. The sign 
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of this term within the troposphere is then determined by 
the increased magnitude and extent of the negative tendency 
due to vertical advection. 

The strong lower tropospheric negative tempera- 
ture tendency, which is present throughout the development of 
this cyclone, is primarily a manifestation of the negative 
advection associated with the upward vertical velocities 
in the storm volume. This temperature decrease is consistent 
with the advective processes affecting the storm volume in 
the lower layers, but the penetration of the large negative 
temperature tendency to 200 mb after 1318 would cause exces- 
sive local cooling of the column in the absence of diabatic 
processes. Excessive cooling of the column is inconsistent* 
with continued cyclogenesis, as mentioned previously. The 
magnitude of this negative upper- layer feature has probably 
been overestimated by the use of the kinematic omega in the 
budget calculation. 

The advective time section for the outer volume 
(Fig. 30b) shows significant reduction in the rate of tem- 
perature increase in the upper layers above 400 mb. This 
change in the heating pattern between the volumes suggests 
that a greater amount of warm advection is concentrated in 
the upper layers at the inner budget volume. 

c. Thermodynamic Energy Budget Diabatic Residual 

The residual temperature change time section for 
the inner budget volume (Fig. 31a) contains several inter- 
esting features. Strong lower tropospheric diabatic effects 
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are indicated at the first two time periods (1206 and 1218) , 
which is consistent with the initial, shallow cyclone growth. 
This warming feature attains its maximum value at 775 mb 
during the second time period. The distribution of this 
heating is consistent with the strong sensible heating and 
shallow convection during the formative period of this cyclone. 

A pause in residual diabatic warming occurs during 
the third time period (1306) . A small amount of cooling is 
indicated in a layer extending from approximately 500-250 
mb. Cooling at these levels during this early stage of cy- 
clone growth is physically unreasonable. A similar pause 
was noted in the vorticity tendencies by Calland (1983), 
which suggests that a problem exists at this time period in 
the ECMWF FGGE analyses. 

Diabatic heating is re-established in the column 
by 1318, except in the lowest layer and above 250 mb. The 
magnitude of this warming is less than during the early 
stages of development, but the vertical growth of the feature 
properly reflects the deeper convection which is present during 
this period. Although the sign of diabatic heating is cor- 
rect, the magnitude appears to have been underestimated, 
considering the amount of convective and large-scale activity 
associated with the system at this time. 

The diabatic cooling feature located in the lower 
layers after 1306 is not physically reasonable. The budget 
determination of this cooling results, in part, from inaccu- 
racies in the estimation of the lower tropospheric advection 
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tendency as the storm transits into a colder air mass between 
1318 and 1406. Although thermal gradients are properly 
resolved at the evaluation times, the 12-h time average of 
the advection term probably underestimates the true hori- 
zontal advection during this period. This artificial decrease 
in advective cooling after 1318 leads to a less positive 
diabatic residual in the budget formulation. The 12-h time 
resolution of the data is not sufficient to resolve the strong 
thermal advective effects when the cyclone moves rapidly into 
the cold air. 

During the last time period (1418) , strong diabatic 
effects with a maximum at 600 mb are indicated in Fig. 31a. 

The restriction of this warming feature to the mid-troposphere 
is consistent with the diabatic process as the cyclone peaks 
in intensity and then occludes. The magnitude of this feature 
(35 K/day) seems to be excessive, but may be a reflection of 
the intensity reached by the cyclone during this period. 

Outer budget volume results (Fig. 31b) are again 
similar to those of the inner volume, except for the smoothing 
associated with the larger area. 

d. Column-Averaged Results 

The area-averaged, column-averaged temperature 
changes due to the diabatic residual and advective tendency 
terms of the thermodynamic energy budget for the inner and 
outer storm volumes are compared in Fig. 32. During the period 
00 GMT 12 January 1979-00 GMT 13 January 1979, the storm 
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center moves primarily to the east and is relatively near 
the Asian east coast conventional, rawinsonde and other data 
sources. The estimation of the horizontal advection tendency 
is relatively accurate as there is little thermal gradient 
in the zonal direction, therefore the diabatic residual 
tendency calculation is believed to be most accurate during 
the first two time periods. Residual diabatic effects reach 
a maximum of 8.2 K/day in the column at storm radius 6, and 
7.6 K/day at storm radius 10 by 1218. This warming maximum 
exceeds the maximum advective cooling (-7.7 K/day at radius 
6 and -6.8 at radius 10) which is not registered until 1406. 

The pause in diabatic warming noted in the verti- 
cal time sections at 1306 is reflected in both the diabatic 
residual and adiabatic results for this time. Residual dia- 
batic effects suffer a sharp, physically unrepresentative 
decrease between 1218 and 1306. This large decrease in magni- 
tude is most probably the result of a problem in the ECMWF 
analyses during this period as previously mentioned. 

One important inference can be drawn from the 
results of the earlier, relatively uncontaminated time periods. 
These results indicate that the magnitude of the contribution 
of diabatic heating to the total column temperature change 
is as large or larger than the advective term, which suggests 
that the diabatic effects are as important as the combined 
effects of horizontal and vertical advection in maritime 
cyclogenesis . 
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C. DIABATIC EFFECTS ESTIMATED BY NOGAPS 



Independent estimates of the total diabatic heating and 
the component terms generated by the NOGAPS diagnostics 
package are interpolated to the storm volume and area-averaged. 
Vertical time sections of the weighted temperature changes 
are constructed for each of the terms. The weighted NOGAPS 
diabatic layer temperature changes are column-averaged to 
produce total column temperature tendencies. These tenden- 
cies are compared directly with the thermodynamic energy 
budget results. 

1 . Procedure 

Diagnostic fields of layer temperature change from 
total diabatic heading and the component terms of convective 
condensation heating, large-scale condensation heating, eddy 
transport of sensible heating, long-wave radiation heating 
and short-wave radiation heating are generated by NOGAPS 
from the ECMWF FGGE analyses. The temperature tendencies at 
the six sigma levels in the NOGAPS model are assumed to be 
representative of heating in the layers centered on these 
levels. The sigma value at each interpolated gridpoint is 
converted to a pressure value using the relationship 

P = o (P - 50) + 50 
s 

where 50 mb is the top of the atmosphere in the NOGAPS model, 

P is the pressure at a given sigma level, and P s is the sur- 
face pressure. Pressure layers are constructed around the 
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converted sigma levels at each gridpoint, and weighting fac- 
tors are computed to normalize these layers to 100 mb, as 
before. The pressure layers and associated weighting factors 
from time 1312 are listed in Table II for illustration. 

These weighting factors are applied to the interpolated 
NOGAPS gridpoint temperature change values at their respec- 
tive levels, and the result is area-averaged over the storm 
vol ume . 

The NOGAPS temperature changes, which are generated 
at 00 and 12 GMT, must be time-averaged to be consistent with 
the thermodynamic energy budget (06 and 18 GMT) . Vertical 
time sections are constructed for -the total diabatic tempera- 
ture tendency and each of the component heating terms. These 
time sections indicate the magnitude of diabatic heating 
that is applied to the storm volume by layer during cyclone 
development. Column averages are calculated from the layer- 

weighted temperatures as discussed in the previous section. 

2 

Total diabatic input (W/m ) to the storm volume is 
calculated for the NOGAPS results using the column-averaged 
diabatic temperature tendencies at radii 6 and 10 for 06 and 
18 GMT each day. The energy input from NOGAPS latent heat 
release in the storm volume is calculated separately for com- 
parison with literature. The energy input to the storm 
volume is determined by 

AT 1 o • o 

(K/day) x c p (J/Kg-K) x 8 ~ g 4 QQ ( da Y/s) * mass (Kg/m ) = q (W/m ) 
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where AT/At is the column-averaged temperature change due to 
the diabatic heating, and the mass is determined for the 
storm volume from the area-averaged surface pressures. This 
mass is adjusted for the model atmosphere which has its top 
at 50 mb. The resulting energy flux values are discussed in 
the next section. 

2 . Analysis 

a. Diabatic Time Sections 

The area-averaged, weighted NOGAPS total diabatic 
fields for storm radii 6 and 10 are plotted as time sections 
in Fig. 33. The overall warming feature for radius 6 appears 
to present a physically reasonable depiction of the diabatic 
process as it occurs over the growth of this intense cyclone. 
During the early stage of cyclone development, warming is 
confined to the lower layers', with cooling indicated above. 

A similar pattern is present in the diabatic residual from 
the temperature budget analysis (Fig. 31a) at times 1206 and 
1218. The warming shown in Fig. 33a increases in magnitude 
and vertical extent with time, until a maximum in both is 
reached at 1318 after the explosive development has been 
established. The rapid upward extension of diabatic warming 
through four sigma levels between 1218 and 1306 is consistent 
with the intense, 12-h vertical development which organizes 
around the storm center just prior to the onset of explosive 
deepening. The vertical retreat of this warming feature 
between 1406 and 1418 indicates that the model diagnoses a 
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weakening in the diabatic process during this period; and 
a termination by 00 GMT 15 January 1979. The maximum of 
diabatic warming leads the peak in storm intensity, which 
occurs at approximately 1800 GMT 14 January 1979, by 18 to 
30 hours. This phase difference between the diabatic heating 
and cyclone intensification appears to be reasonable if the 
diabatic forcing is responsible for, or at least is contribut- 
ing to, the atmospheric response. 

The diabatic effect that is diagnosed for the 
outer storm volume in Fig. 33b is very similar to that of the 
inner storm volume, except that it has been smoothed and re- 
duced in magnitude from averaging over a larger region. This 
similarity in the diabatic term between the storm volumes 
indicates that the most significant diabatic heating is not 
confined to the inner storm volume, but is spread over a 
much larger region around the storm center. The only signi- 
ficant difference in the heating time sections between the 
storm volumes is in the location of the warming maximum at 
radius 10, which has been shifted forward by 6 hours to time 
period 1318. This shift is later shown to be the result of 
the distribution of the sensible heat flux. 

Individual contributions to the diabatic term 
are evaluated in the time section format. The vertical con- 
tribution of model-diagnosed convective precipitation, 
large-scale precipitation, eddy transport of sensible heating 
and long-wave radiation to total diabatic heating are 
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included in Figs. 34-37. Short-wave radiation time sections 
are not shown in the figures, but are discussed below. 

The convective condensation heating field includes 
both deep and mid-level convection. After time period 1218, 
this heating input in the inner volume (Fig. 34a) determines 
the bulk of the total diabatic temperature tendency present 
in Fig. 33a. The convective contribution to the outer volume 
is similar, except that more heating is indicated during the 
first two time periods. This early increase in heating is 
due to unrelated convection to the south of the storm center 
which is included in the larger budget volume. 

The model properly attempts to center the maxi- 
mum of convective heating in the lower troposphere during 
time period 1318. It establishes this maximum just above 
sigma level 6 (792 mb) at both storm radii, which is not 
unreasonable, but may be somewhat low considering the deep 
convection which develops in the vicinity of the storm center 
during this period. 

The vertical eddy transport term makes its most 
significant contribution to total diabatic heating in the 
lowest sigma layer within storm radius 6, as illustrated in 
Fig. 35. The bulk of the diabatic heating at the inner 
volume during the first two time periods is determined by 
this term. This early influence of the sensible heat flux 
on total diabatic heating within the inner storm volume is 
not unexpected, since the storm center is embedded in a 
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strong surface sensible heat flux feature from 12 GMT 12 
January 1979 to 00 GMT 13 January 1979. 

The cooling feature present in the lowest layer 
during the last two time periods in Fig. 35a results from a 
model-diagnosed, downward sensible heat flux in advance of 
the storm center. As described in the previous chapter, the 
magnitude of this feature is physically unreasonable. In 
combination with the long-wave radiational effects, this 
feature accounts for the inaccurate diagnosis of low-level 
diabatic cooling at radius 6 during this period (Fig. 33a) . 

The vertical eddy transport time section for the 
outer storm volume (Fig. 35b) is very similar to that of 
the inner volume. The addition of sensible heat to the storm 
volume is extended in time to 1406, as the larger area covered 
by this storm radius encompasses more of the sensible heat 
flux features to the west of the storm center. The effect 
of the sensible cooling during the last two time periods is 
greatly reduced within radius 10, as the sensible cooling 
maximum is primarily confined to the inner budget volume. 

Large-scale condensation heating that is diagnosed 
by NOGAPS (Fig. 36) contributes little to the total diabatic 
term. The relative insignificance of large-scale heating 
initially seems unreasonable, when it is considered that this 
cloud source accounts for the bulk of the extensive cloud 
mass which is located to the northeast of the storm center 
throughout most of the evaluation sequence. However, the 
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magnitude of the vertical motion is small until explosive 
deepening is established at 1318. By this time,, the storm is 
in relatively cool air where the saturation vapor pressure 
is low, which restricts the magnitude of the potential con- 
densation heating. 

The vertical time section from the inner storm 
volume (Fig. 36a) indicates that the effect of large-scale 
heating is confined to the lowest sigma layer until explosive 
cyclogenesis is established at 1318. Large-scale condensation 
heating then extends through four sigma layers, which is 
consistent with the diagnosed deepening of the clouds to the 
northeast of the storm center during the explosive stage of 
development. However, the temperature tendencies associated 
with this term are an order of magnitude less than the con- 
vective term, which illustrates why deep convection is so 
important in cyclogenesis. The outer storm volume (Fig. 36b) 
provides similar results, although the heating maximum is 
shifted forward by 12 h to 1406. The upper level contribution 
is smoothed out as it is averaged over the larger storm area, 
but the contribution in the lowest layer is increased slightly 
at the later time periods as the more extensive low-level 
cloudiness is included in the larger storm volume. 

The long-wave radiation time sections (Fig. 37) 
reflect cloud-top cooling and atmospheric IR cooling. The 
most significant long-wave cooling occurs in the lowest sigma 
layer where the most extensive cloudiness is located. This 
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cooling is strongest during the last two time periods, and 
provides a weighted contribution of -9 K/day to the low-level 
cooling feature in the total diabatic temperature tendency 
time section at radius 6 (Fig. 33a) . The strong, low-level, 
long-wave radiative cooling contribution as determined in 
the NOGAPS radiation parameterization infers a cloud-top 
instability mechanism similar to that described by Petterssen, 
et al . (1962). Vertical time sections of short-wave radiation 

heating are not shown. Diagnosed heating contributions from 
short-wave radiation are negligible (weighted maximum of 0.2 
K/day) and have little effect on the total diabatic term, 
b. Column Averaged Results 

The column-averaged, storm volume temperature 
tendencies are computed for total diabatic heating and the 
component terms. These results are compared for the inner 
and outer storm volume in Fig. 38. As noted in the evaluation 
of the NOGAPS time sections, a remarkable similarity exists 
in the display of the results for radius 6 and radius 10, 
which reflects the large storm volume affected by similar 
diabatic processes. 

A nearly linear increase in the total diabatic 
term is indicated in Fig. 38 until a maximum is reached when 
explosive deepening is established at 1318. A linear decrease 
in this warming then follows. This representation of diabatic 
heating is very reasonable. The maximum diabatic column 
temperature change rate of 6.96 K/day at radius 6 compares 
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favorably with results from Yanai, et al. (1973) and Lord 
(1982) who showed average column temperature changes from 
diabatic heating in tropical cloud clusters of approximately 
5 K/ day . 

An inspection of Fig. 38 reveals that latent heat 
release accounts for almost all of the net diagnosed warming. 
Only a small portion of this latent heat release is attrib- 
uted to large-scale cloud processes. The vertical eddy 
transport of sensible heat from the surface makes a signifi- 
cant contribution to diabatic warming (radius 6) during the 
early time periods, but reverses in sign at time period 1406 
when (questionable) sensible cooling is diagnosed to the 
southeast of the storm center. Much of the column warming by 
the surface flux during the early period of cyclone develop- 
ment is offset by long-wave radiative effects, which cool 
the column by 1. 5-2.0 K/day . Short-wave radiative warming 
of the column is less than 0.02 K/day and is not plotted in 
the figure. 

c. Energy Input from NOGAPS Latent Heat Release 

The diabatic energy input to the storm volume 
from NOGAPS-diagnosed latent heat release is computed for the 
inner and outer storm radii. A comparison of the energy 
input within the inner and outer storm volumes (Fig. 39) 
indicates that the level of input remains relatively un- 
changed between the volumes. Differences at the first time 
period are attributed to the inclusion of an unrelated cloud 
mass to the south in the outer storm volume. 
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The diabatic energy input from latent heat re- 

2 

lease for radius 6 has a maximum of 841 W/m at time 1318, 

2 

with an average input for all time periods of 352 W/m . These 

values differ radically with those determined from case 

studies of continental cyclogenesis by Danard (1966) and Bullock 

and Johnson (1971) , who found the average storm volume energy 

2 2 

input from latent heat release to be 11 W/m and 8 W/m , 

respectively. A much more favorable comparison is made with 

results from Yanai, et al . (1973) for latent heating in 

tropical cumulus cloud clusters over the Marshall Islands. 

2 

They determined an average energy input of 485 W/m from 
latent heating in the column. This comparison underscores 
the enhanced role of diabatic heating in maritime cyclogene- 
sis. In relation to latent heating, the energy input from 

net diabatic heating for radius 6 reaches a maximum of only 
2 

760 W/m at 1318. The average net diabatic energy input to 

2 

the inner storm volume for all time periods is 266 W/m . 

D. NOGAPS AND BUDGET COLUMN WARMING COMPARISON 

A comparison of NOGAPS total diabatic warming and the 
thermodynamic energy budget residual diabatic warming 
(Fig. 40) indicates that the magnitudes of the maxima are 
similar (8.2 K/day vs. 6.96 K/day at radius 6) , but the NOGAPS 
maximum lags that from the budget diagnosis by 24 h. This 
comparison is not completely fair since problems are indicated 
with the budget calculation, and the budget residual 



91 



diabatic warming probably did not peak during time period 
1218, as shown. In any event, the temperature tendency 
maxima appear to be of the same magnitude, and the timing 
of the NOGAPS diabatic heating seems to be more physically 
representative. A similar comparison holds for the outer 
storm volume (Fig. 40b) . It seems that the NOGAPS diabatic 
package responds to the large-scale forcing of the FGGE data 
much better than the thermodynamic energy budget over the 
open ocean. This difference in the quality of heating esti- 
mates is probably connected to the sophisticated BL physics 
included in the NOGAPS diagnosis, but missing from the budget 
analysis. The NOGAPS package permits a direct evaluation 
of the diabatic heating over the ocean while a budget compu- 
tation can only infer this process indirectly. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The contribution of diabatic heating in a case of ex- 
plosive cyclogenesis in the western North Pacific Ocean during 
FGGE was investigated in this thesis. An evaluation of the 
NOGAPS heating package was conducted using the ECMWF analyses 
for this cyclone during 13-15 January 1979. The storm-related 
thermal budget was computed with the application of quasi- 
Lagrangian diagnostics to the analyses. The diabatic tem- 
perature tendencies from the NOGAPS diagnostics and the 
temperature budget were examined and discussed. 

A. RESULTS 

1 . NQGAPS Boundary Layer and. Cloud Parameterization 
The model boundary layer parameterization yields 
correct positioning of the sensible heat flux features 
relative to computed air-sea temperature differences and the 
boundary layer (ECMWF 1000 mb) wind fields. Generally, 
positive flux features are confined to the cold air behind 
the front and negative fluxes are located in the warm sector 
to the east of the front. The large negative flux maxima 
generated by the model at later time periods are a reflection 
of the ECMWF wind and temperature analyses and not the NOGAPS 
model performance. The magnitudes of the diagnosed positive 
sensible heat flux maxima are high, but not excessive com- 
pared to previous studies given the intense nature of this 
cyclone . 



The top of the model planetary boundary layer is 
collapsed to extremely low heights in regions where sensi- 
ble heat flux maxima are diagnosed. Modification of the drag 
coefficients in these regions may lead to an artificial 
increase in the sensible heat flux values in the model. 

The combined cloudiness diagnosed in the model is 
representative of the verifying satellite imagery. Outbreaks 
of convection observed in the vicinity of the storm center, 
as well as observed deep frontal convection, are handled well 
by the model parameterization. A tendency to displace these 
features slightly to the west of their verifying positions 
is noted at several of the evaluation times. Cloud streets 
and open-cell cumulus clouds in the cold air streams are 
treated as unstable stratus in the parameterization. This 
model-generated, large-scale cloud feature is. diagnosed over 
regions of strong sensible heat flux. The boundary layer is 
collapsed at gridpoints where these clouds are diagnosed in 
the boundary layer. The model depiction of the large-scale 
cloudiness properly accounts for the large cloud mass to the 
northeast of the storm center, as well as a portion of the 
frontal clouds, throughout most of the storm sequence. 

Model-diagnosed temperature tendencies associated with 
large-scale condensation heating are reasonable. Temperature 
tendencies associated with convective condensation heating 
appear to be excessive, but are related to the very large 
convective precipitation rates determined by the model for 
this system. These precipitation rates are not totally 
inconsistent with the observed intensity of the cyclone. 
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2 . NOGAPS and Thermodynamic Energy Budget Diabatic Heating 

The maximum of NOGAPS-diagnosed diabatic heating 
precedes the peak in storm intensity by approximately 24 h. 

This temporal distribution of the heating is a reasonable 
reflection of the explosive development of the cyclone and 
the time lag required between atmospheric forcing and atmos- 
pheric response. 

NOGAPS-diagnosed total diabatic heating is primarily 
composed of latent heat release in the column. The great 
bulk of the latent heating is convective in nature, with 
large-scale processes accounting for only a small percentage 
of the heating. All other NOGAPS heating terms are small in 
comparison to the latent heat release. 

The NOGAPS average diabatic energy input to the storm . 
volume through latent heat release is an order of magnitude 
larger than estimates based on continental cyclogenesis. The 
results are in general agreement with a previous study of 
latent heating in maritime, tropical cumulus cloud clusters. 

The large, diagnosed contribution of latent heat release to 
the energy of the system in this case of maritime cyclogene- 
sis underscores the increased importance of the diabatic 
processes in cyclogenesis over the oceans. 

A significant amount of diabatic heating is determined 
from thermodynamic energy budget calculations during the first 
two time periods (1206-1218), as the cyclone experiences 
initial growth. The magnitude of the budget diabatic residual 
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decreases unrealistically as the storm moves away from land 

after the third time period (1306) . This calculated deteriora- 

i 

tion of the diabatic process coincies with the explosive 
development of the system, when the contribution from the 
diabatic term is expected to be greatest. The magnitude of 
diabatic heating increases significantly during the last time 
period (1406) when the system reaches peak intensity and 
the storm center is again near land where the quality of the 
data is much better. 

Errors occur in the thermodynamic energy budget where 
the storm tracks over the open ocean. The temperature time 
tendency and advective terms of the thermodynamic energy 
equation do not appear to be properly represented using the 
12-h, large-scale FGGE forcing as the storm center moves 
away from land. 

Adiabatic temperature advection produces net cooling 
in the column at all times. The magnitude of the diabatic 
warming is larger than this cooling during the early stages 
of cyclone development when the diabatic calculation appears 
to be relatively unaffected by budget error. 

Differences in the magnitude of the maximum diabatic 
column temperature change as determined by NOGAPS and the 
thermodynamic energy budget are less than 15%, with the 
thermodynamic energy budget producing the greater maximum. 

The apparent 24 -h lag between the maxima appears to be arti- 
ficially induced by discrepancies in the ECMWF analyses and 
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the smoothing effect of the 12-h time averaging of the 
advective terms in the budget formulation. 

Maximum diabatic forcing in NOGAPS precedes the 
maximum advective cooling calculated in the thermal budget 
by 12 h. The comparison indicates that the NOGAPS-diagnosed 
diabatic heating is of the same magnitude as, and maintains 
the proper temporal distribution with, the calculated advec- 
tive cooling. The similarity of results from the inner and 
outer storm volumes indicate that the diabatic effects are 
spread evenly over a relatively large region around the storm 
center . 

The results of this study indicate that the large-scale 
forcing by the FGGE data is sufficient to produce reasonable 
heating estimates from the NOGAPS diabatic package, but is 
not sufficient to adequately resolve the diabatic heating 
through application of a thermodynamic energy budget calcu- 
lation over the open ocean. 

B. RECOMMENDATIONS FOR FUTURE RESEARCH 

Further investigation of NOGAPS boundary layer behavior 
in regions of strong sensible heating is required to evalu- 
ate the effect of the collapsing boundary layer and resulting 
modification of the drag coefficients on the diagnosis of 
the surface sensible heat flux. 

The thermodynamic energy budget calculation and NOGAPS 
diagnostics should be applied together in other cyclone 
evaluations to verify the apparent superiority of the NOGAPS 
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diabatic heating estimates over open ocean regions which was 
observed in this study. The flux form of the thermodynamic 
energy equation should be applied in the budget to see if any 
improvement in the calculations result. 

The thermodynamic energy budget calculation and NOGAPS 
diagnostics should be performed using temperatures derived 
from the ECMWF height fields to determine if the initialized 
temperatures were a large source of error in the calculations . 
A storm-area moisture budget should be determined as the next 
phase in the evaluation of this storm system. 

C. FINAL CONCLUSION 

Diabatic heating produces significant storm-area column 
temperature changes in explosive maritime cyclogenesis, which 
enhances system intensification. Proper model parameteriza- 
tionsof these processes are essential if accurate numerical 
weather forecasts are to be realized over the oceans. 
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APPENDIX A 



SYNOPTIC SYMMARY 

A complete synoptic discussion of the western Pacific 
Ocean storm system evaluated in this thesis is given by 
Calland (1983) . The significant events of cyclone develop- 
ment are summarized by time period to present an overview of 
cyclogenesis as it occurred during 12-15 January 1979. 

The large-scale flow pattern at 00 GMT 12 January 1979 
transports cold, continental air southward into the vicinity 
of a warm, southerly flow from the tropics. A reg-ion of 
strong, low-level baroclinity develops to the east of Japan. 

A quasi-stationary , long-wave trough is located over the 
region in the mid-troposphere. This trough supports a steady 
flow of cold air off the Siberian coast. The mid-level flow 
over the area of the incipient system is primarily zonal. A 
broad jet is present in the upper troposphere, and reflects 
the baroclinity in the low levels. A jet streak (> 15 m/s) 
remains well upstream of the region of the developing 
disturbance . 

A. 00 GMT 12 JANUARY 1979-00 (24T 13 JANUARY 1979 

The incipient storm has a closed center of 995 mb. A 
weak thermal ridge is established near the surface center 
and indicates continued f rontogenesis . The combination of 
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flow from a strengthening Siberian High to the north and a 
nearly stationary cold low to the east pumps cold air 
southward over the warmer sea surface, which continues the 
destabilization of the surface layer. A polar low has formed 
to the west of the incipient system over the southeast tip 
of Japan. 

In the mid-troposphere, a short wave which incorporates 
a broad area of positive absolute vorticity has formed about 
800 n mi upstream of the incipient system. A small region of 
positive vorticity has formed over the storm center. The jet 
core in the upper level is superposed over the surface low 
center and the jet streak continues to approach from the 
west . 

B. 00 GMT 13 JANUARY 1979-00 GMT 14 JANUARY 1979 

The storm enters the explosive stage during this period, 
as the central pressure falls to 977 mb. The surface low 
accelerates to the northeast into the vicinity of the maximum 
sea-surface temperature gradient. A diffluent trough begins 
to form in the mid-troposphere. Mid-level vorticity increases 
as the jet streak propagates into the region from the west 
and enters the diffluent trough. 

The storm system experiences increased vertical develop- 
ment following the jet streak interaction. The arrival of 
the jet streak at 250 mb (max > 85 m/s) over the surface 
system enhances low-level circulation. The diffluence and 
cyclonic shear that are associated with this streak increase 
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both vertical motion through the mid-troposphere, and storm 
circulation in general. The short-wave that was present over 
Japan during the previous period is now absent, indicating 
that further explosive development is not sustained by short- 
wave, positive vorticity advection aloft. 

C. 00-12 GMT 14 JANUARY 1979 

The central pressure falls to 965 mb during this 12-h 
period. The amplitude of the surface thermal ridge increases, 
indicating enhanced circulation in the lower troposphere. 
Self-amplification continues at 500 mb. The apparent lack 
of significant vorticity advection indicates that the con- 
tinued intensification of the cyclone system is not forced 
in the mid-troposphere. Cyclonic shear and curvature to the 
left of the jet axis increases cycloni-c vorticity at 250 mb 
and enhances low-level storm circulation. Vertical velocities 
increase and cyclone development extends to higher levels 
in the troposphere. 

D. 12 GMT 14 JANUARY 1979-00 GMT 15 JANUARY 1979 

The first half of this period marks the most explosive 
stage of cyclone development. The surface pressure falls an 
additional 18 mb and the cold air begins to wrap around the 
vortex, which signals the onset of occlusion. The mid-level 
trough continues to amplify, but its vorticity center elon- 
gates. This change in the vorticity pattern foreshadows 
weakening of the system. The jet stream is still in a 
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favorable position over the surface low to promote continued 
intensification, but the rapid movement of the cyclone to 
the northeast will reduce this jet influence. 

E. 00 GMT 15 JANUARY 1979 

The surface low is located over the Aleutian Islands at 
this time. The central pressure reaches 947 mb as the system 
begins to occlude. A closed circulation has formed at 500 mb 
and negative vorticity advection into the region induces 
subsidence and filling. The jet has decreased in intensity 
with the weakened low-level baroclinity, and the cyclone 
center is located to the north of the jet axis. 
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APPENDIX B 



FGGE DATA 

The data base used in this study is the First GARP Global 
Experiment (FGGE) Level Ill-b objectively analyzed meteorologi- 
cal fields. These analyses were produced by the European 
Center for Medium-range Weather Forecasting (ECMWF) from the 
first Special Observing Period (SOP-1) for the region cover- 
ing the track of the western Pacific cyclone of interest. 

This data set is the product of an intensive one year world- 
wide data collection effort which utilized conventional sur- 
face, aircraft, rawinsonde , pilot balloon, drifting buoy and 
satellite observations to produce global meteorological fields 
for research purposes. The data collection for this project 
was carried out under the auspices of the Global Atmospheric 
Research Program (GARP) for the year 1979. SOP-1 covered 
the first two months of the project (5 January — 5 March 1979) 
and is significant in that more observation platforms were 
brought to bear during this period with the resulting analyzed 
fields archived in 6 vice 12-h intervals. 

In FGGE level Ill-b, the processed temperature, moisture, 
wind velocity and height data have been transformed into 
dynamically-derived meteorological fields through a four- 
dimensional data assimilation technique described by 
Bengtsson et al . (19 82) . These fields were produced in 15 

vertical levels with a horizontal resolution of 1.875° in 
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latitude and longitude. The data was error checked at the 
ECMWF by comparison of first guess fields and nearby observa- 
tions level by level. As reported by Bengtsson et al . (1982) 

the SOP-1 data was generally found to be of excellent quality 
The fields available in the FGGE data are compiled by level 
in Table III. 

Validation of the ECMWF analyses has been performed in 
several studies including that of Calland (1983) in his 
assessment of the analyses used in this thesis. Through a 
comparison of the ECMWF level III-b, NMC and FNOC analyses 
during the period of storm development, Calland found the 
ECMWF analyses to be very good overall. He did note, however 
that the 0600 and 1800 GMT analyses contained spurious poten- 
tial temperature biases which adversely affected his budget 
computations. 

FGGE vertical velocity and temperature fields were taken 
from the initialization stage of the data assimilation at 
the ECMWF, rather than from the analysis stage. The magni- 
tudes of the resulting FGGE vertical velocities were small 
compared to those determined kinematically. Kinematic 
vertical velocities are substituted for the FGGE values in 
the data set utilized in this evaluation. The FGGE tempera- 
tures are verified against the height (mass) field to ensure 
consistency . 
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APPENDIX C 



THE NAVY OPERATIONAL GLOBAL ATMOSPHERIC PREDICTION SYSTEM 

(from Ranelli, 1984) 

The Navy Operational Global Atmospheric Prediction Sys- 
tem (NOGAPS) used at the Fleet Numerical Oceanography Center 
(FNOC) is a slightly modified version of the UCLA general 
circulation model. NOGAPS has been the Navy's operational 
atmospheric forecast model since August 1982. The following 
sections describe the various features of NOGAPS as used 
during the experiment. The complete model has been des- 
cribed by Rosmond (1981) . 

A. DYNAMICS 

The dynamics of the UCLA GCM are described in detail by 
Arakawa and Lamb (1977) and are only discussed briefly here. 
NOGAPS is a primitive equation model. The- prognostic varia- 
bles are horizontal velocity, V, temperature, T, surface 
pressure, p^ , and specific humidity, q. Additional prognos- 
tic variables associated with the planetary boundary layer 
(PBL) will be described below. The finite difference scheme 
used has a spatial resolution of 2.4° lat by 3.0° long. The 
variables are staggered in the horizontal according to Arakawa 
scheme C. The center grid point contains the T value. The 
meridional wind component, v, is carried at points north and 
south of the center point and the zonal wind component, u, is 
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carried at points east and west of the center point. The 
numerical differencing scheme is both energy and enstrophy 
conserving . 

NOGAPS uses a sigma coordinate system in the vertical 
defined as: 



a = (p - p i ) /it 



where : 



p. = 50 mb and tt = p - p. , 

p is pressure and p g is surface pressure. There are six 
model layers in the vertical with the top of the model atmos- 
phere at 50 mb. All prognostic variables except vertical 
velocity, a, are carried at the middle of each layer - . Verti- 
cal velocity is carried at the layer interfaces. 

NOGAPS uses a second order (leapfrog) time difference 
scheme with a four minute time step. Model diabatics are 
executed every forty minutes. A Matsuno time step is used 
every fifth time step. This is used to control the compu- 
tational mode and to assist in the assimilation of the diabatic 
effects. In regions above 60° latitude, a special Fourier 
filter is used to avoid an extremely short time step. Whereas 
a simple three point filter is used equatorward of 60 deg. 

This filtering reduces the amplitudes of the zonal mass flux 
and pressure gradients and maintains computational stability. 



B. MODEL DIABATICS 



The sophisticated model diabatics contained in NOGAPS is 

an important component in this experiment. This treatment 

of the diabatic processes is necessary to adequately simulate 

fluxes across the air-sea interface and to propagate the full 

effect of these changes throughout the atmosphere. NOGAPS 

directly computes the physical processes for: 

dry convective adjustment 
large-scale precipitation 
diagnosis of stratus cloud depth 
mid-level convection 
ground hydrology 
surface friction 

horizontal diffusion of momentum 
radiative transfer processes 
cumulus convection 

1 . Planetary Boundary Layer 

The planetary boundary layer (PBL) is defined as a 
well mixed layer in moisture, moist static energy and momen- 
tum. It is assumed to be capped by discontinuities in tem- 
perature, moisture and momentum. The PBL treatment in this 
model follows Deardorff (1972) and has been formulated for the 
UCLA GCM by Randall (1976) . It allows for interaction be- 
tween the PBL and cumulus cloud ensembles and/or a stratus 
cloud layer at each grid point. Surface fluxes are determined 
using a bulk Richardson number based on the values of the sea 
surface temperature and the values of V, T and q from the 
adiabatic portion of the model. These values are then used 
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to predict a new PBL depth and the strength of the inversion 
j umps . 

The NOGAPS PBL is constrained to remain in the bottom 
sigma level of the model. This differs from the original 
formulation of the UCLA GCM, in which the PBL was allowed 
to pass out of this layer. An overly deep PBL can result 
in serious computational problems with the model. Constrain- 
ing the PBL this way imposes a maximum depth of about 200 mb 
on the PBL. 

2 . Cumulus Parameterization 

Cumulus parameterization in NOGAPS follows the scheme 
of Arakaw a- Schubert (1974) as introduced into the model by 
Lord (19 78) . In the model, cumulus clouds must have the^ir 
bases at the top of the PBL. Cloud tops can be at all sigma 
levels above the PBL. Cumulus clouds are modeled as entrain- 
ing plumes in which environmental air is mixed with the PBL 
air from which the cloud originated. Tendencies of moisture, 
temperature and momentum are diagnosed as well as the cloud 
mass flux. The cloud base mass flux removes mass from the 
PBL, which decreases the PBL depth. Condensation occurs at 
each grid point where the air becomes supersaturated. A 
moist convective adjustment procedure removes convective 
instability between mid-tropospheric layers that is not 
eliminated by clouds originating from the PBL. 

3 . Radiation 

The radiation parameterization follows Katayama 
(1972) and Schlesinger (1976). It includes both a diurnal 
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variation and interaction with the cloud distribution. 

Radiative transfer processes for incoming solar radiation 
are computed. Effects of water vapor, Rayleigh scattering 
by air molecules and absorption and scattering by water 
droplets in clouds are included. Reflection due to clouds is 
also calculated. The model cloud cover predicted by the PBL, 
the cumulus parameterization and large-scale precipitation 
interact with the long wave radiation. The net surface heat 
flux is computed as a function of the incoming solar heat 
flux, long-wave radiation and sensible heat flux. In the 
present model, this affects only the surface temperature over 
bare land and ice and has no effect on sea-surface temperature. 
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APPENDIX D 



TABLES 



TABLE I 

SENSIBLE HEAT FLUX COMPARISON 



INVESTIGATOR 


MAX 

FLUX 

W/m^ 


EVALUATION 

TIME 


MAX 

FLUX 

W/m^ 


Petterssen 

(1962) 


699 


13/0000 


1400 


Pyke (1964) 


105 


13/1200 


900 


Manabe (1958) 


498 


14/0000 


1100 


Ga 11 Sc 

Johnson (1970) 


1080 


14/1200 


500 


Bosart (1984) 


800 


15/0000 


300 
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TABLE II 



SIGMA AND PRESSURE LAYER RANGES AND WEIGHTING 
TIME PERIOD 12 GMT 13 JANUARY 1979 



SIGMA LAYER 


PRESSURE LIMITS (MB) 


WEIGHTING FACTOR 


6 


994-696 


3.0 


5 


696-497 


2.0 


4 


497-298 


2.0 


3 


298-149 


1.5 


2 


149-75 


0.74 


1 


75-50 


0.24 


PRESSURE LAYER 


LIMITS (MB) 


WEIGHTING FACTOR 


1 


994-850 


1.5 


2 


850-700 


1.5 


3 


700-500 


2.0 


4 


500-400 


1.0 


5 


400-300 


1.0 


6 


300-250 


0.5 


7 


250-200 


LO 

• 

o 


8 


200-150 


LO 

• 

o 


9 


150-100 


0.5 



111 



TABLE III 



AVAILABILITY OF FGGE DATA 
(From Bengsston, 1982) 







z 


u 


V 


T 


w 


RH 


10 


mb 


X 


X 


X 


X 


X 




20 


mb 


X 


X 


X 


X 


X 




30 


mb 


X 


X 


X 


X 


X 




50 


mb 


X 


X 


X 


X 


X 




70 


mb 


X 


X 


X 


X 


X 




100 


mb 


X 


X 


X 


X 


X 




150 


mb 


X 


X 


X 


X 


X 




200 


mb 


X 


X 


X 


X 


X 




250 


mb 


X 


X 


X 


X 


X 




300 


mb 


X 


X 


X 


X 


X 


X 


400 


mb 


X 


X 


X 


X 


X 


X 


500 


mb 


X 


X 


X 


X 


X 


X 


700 


mb 


X 


X 


X 


X 


X 


X 


850 


mb 


X 


X 


X 


X 


X 


X 


1000 


mb 


X 


X 


X 


X 


X 


X 
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APPENDIX E 



FIGURES 




Figure 1. Storm Track and FNOC Sea-surface Temperature 

Analysis for (A) 00 GMT 12 January 1979, (B) 12 

GMT 16 January 1979. Contour Interval is 2°C. 
Dashed Portion of the Track not Covered by the 
Analysis. Label X Identifies 12 h Storm Center 
Positions . 
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Figure 2. NOGAPS BL Diagnoses 12 GMT 13 January 1979 of (A) 

Sensible Heat Flux and (B) LCI. Flux Contour Inter- 
val is 100 W/m^ . Dashed (Negative) Values Indicate 
Flux into the Ocean. LCI Greater than 1.0 Indicates 
Unstable Stratus. Label X Identifies the Storm 
Center . 
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Figure 3. NOGAPS BL Diagnoses 12 GMT 13 January 1979 of 
(A) BL Thickness < 20 mb and (B) Sigma Level 6 
Large-scale Precipitation. Contours > 1.0 Indicate 
Location of Respective Fields. Lable X Identifies 
the Storm Center. 
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Figure 4. BL Wind and Flux Fields for 00 GMT 13 January 1979 

of (A) FGGE 1000 mb Wind Field and (B) Sensible Heat 
Flux Field. Isotach Contour Interval is 10 m/s. 

Wind Barbs Indicate Knots. Flux Unit and Contour 
Interval as in Figure 2. Label X Identifies the 
Storm Center. 
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Figure 5. Similar to Figure 4 Except for 12 GMT 13 January 1979. 
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Figure 6. Similar to Figure 4 Except for 00 GMT 14 January 1979. 
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Figure 7. Similar to Figure 4 Except for 12 GMT 14 January 
1979 and Southwest Corner of Grid Shifted to 
30 °N, 150 °E . 
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Figure 8. Similar to Figure 4 Except for 00 GMT 15 January- 
1979 and Southwest Corner of Grid Shifted to 
30 °N, 150 °E . 
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Figure 9 . 



DMSP Visual Satellite Imagery for 2341 GMT 12 
January 1979. Label A Indicates Incipient System. 
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Figure 10 . 



DMSP Infrared Satellite Imagery for 2341 GMT 12 
January 1979. Label A Indicates Incipient System. 
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Figure 11. NOGAPS Cloud Diagnoses of (A) Convective and 
(B) Large-scale for 00 GMT 13 January 1979. 
Contours Indicate Number of Sigma Levels Where 
Condensation Heating is Diagnosed. 
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Figure 12. 



Similar to Figure 9 Except for 
January 19 79 . 



1405 GMT 13 
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Figure 13. Similar to Figure 10 Except for 1405 GMT 13 

January 1979. Label A Indicates Convection to the 
West of the Storm Center. 
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Figure 14. Similar to Figure 11 Except for 12 GMT 13 
January 1979. 
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Figure 15. Similar to Figure 9 Except for 0023 GMT 14 

January 1979. Label A Identifies Approaching 
Polar Low. 
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Figure 16. Similar to Figure 10 Except for 0023 GMT 14 
January 1979. Label A Identifies Polar Low. 
Label B Indicates Low-level Frontal Clouds. 
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Figure 17. 



Similar to Figure 11 Except for 00 GMT 14 
January 19 79 . 
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Figure 18. 



Similar to Figure 9 Except for 1205 GMT 14 
January 1979. Label A Identifies Secondary System. 
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Figure 19. Similar to Figure 10 Except for 1205 GMT 14 

January 1979. Label A Identifies Secondary System 
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Figure 20. Similar to Figure 11 Except for 12 GMT 14 

January 1979 and Southwest Corner of Grid Shifted 
to 30 °N , 150 °E . 
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Figure 21. Similar to Figure 11 Except for 00 GMT 15 January 
1979 and Southwest Corner of Grid Shifted to 
30 °N , 150 °E . 
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Figure 22. Condensation Heating Fields For 00 GMT 13 January 
1979 For (A) Large-scale and (B) Convective. 

Units are K/day. Convective Contour Interval is 
10.0. Large-scale Contour Interval is 0.5. 

Label X Identifies the Storm Center. 
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Figure 23. Similar to Figure 22 Except for 12 GMT 13 
January 1979. 
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Figure 24. Similar to Figure 22 Except for 00 GMT 14 
January 1979. 
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Figure 25. Similar to Figure 22 Except for 12 GMT 14 January 
1979. Convective Contour Interval is 5.0. 
Southwest Corner of Grid Shifted to 30 °N, 150 °E. 
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Figure 26. Similar to Figure 22 Except for 00 GMT 15 January 
1979. Convective Contour Interval is 5.0. 
Southwest Corner of Grid Shifted to 30°N, 150°E. 
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Figure 27. Distribution of Budget Parameters. (A) Layer 

Averaged and Level Parameters, (B) Time Differen- 
cing and Time Averaging Applied. H.A. Refers to 
Horizontal Advection. V.A. Refers to Vertical 
Advection . 
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Figure 28. Layer Potential Temperature from ECMWF Initiali- 
zation and Height Analysis for 700-500 mb Layer 
at Radius 6. Analysis Time 1 Refers to 0000 GMT 
12 January 1979. Subsequent Times are Incremented 
by 12 h. 
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Figure 29. Temperature Time Tendency Time Sections for 

(A) Radius 6 and (B) Radius 10. Contour Interval 
is 2.5. Units are K/day. Dashed (Negative) 
Values Indicate Cooling. Period 1206 Refers to 
00-12 GMT 12 January 1979. 
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Figure 30 . Advective Temperature Tendency Time Sections 

For (A) Radius 6 and (B) Radius 10. Units and 
Contours as in Figure 29 . 
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Figure 31. Diabatic Residual Temperature Tendency Time 
Section for (A) Radius 6 and (B) Radius 10. 
Units and Contours as in Figure 29. 
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Figure 32. Budget Residual (Solid) and Advective (Dashed) 
Column-Averaged Temperature Tendency for (A) 

Radius 6 and (B) Radius 10 . Advective Values Are 
Multiplied by -1. Units Are K/Day . Time Period 1 
Refers to 1206. Subsequent Time Periods Are 
Incremented by 12 h. 
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Figure 33. NOGAPS Total Diahatic Tendency Time Sections For 
(A) Radius 6 and (B) Radius 10. Units and Con- 
tours as in Figure 29 . Vertical Coordinate 
Represents Sigma Levels. 
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Figure 34. NOGAPS Convective Condensation Temperature 
Tendency Time Sections For (A) Radius 6 and 
(3) Radius 10. Units and Contours as in Fiaure 
29 . 
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Figure 35. NOGAPS Eddy Transport of Sensible Heating Time 
Sections for (A) Radius 6 and (B^ Radius 10. 
Units and Contours as in Figure 29 . 
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Figure 36. NOGAPS Large-scale Condensation Temperature 
Tendency Time Sections For (A) Radius 6 and 
(B) Radius 10. Units as in Figure 29. 
Contour interval is 0.5. 
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Figure 37. NOGAPS Long-wave Radiation Temperature Tendency 

Time Sections For (A) Radius 6 and (B) Radius 10 . 
Units as in Figure 29. Contour Interval is 0.5. 
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Figure 38. NOGAPS Column-averaged Total Diabatic, Latent 
Heating, Eddy Sensible Heating and Long-wave 
Radiation Temperature Tendencies For (A) Radius 
6 and (B) Radius 10. Times and Units as in 
Figure 36 . 
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Figure 39 . 



NOGAPS Latent Heating Energy Input 
and Radius 10. Units are in W/ra^ . 
in Figure 36 . 



For Radius 6 
Times as 
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A 



B 





Figure 40. Column-averaged Budget Residual and NOGAPS 

Total Diabatic Temperature Tendencies For (A) 
Radius 6 and (B) Radius 10. Times and Units as 
in Figure 36. 
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